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Notation 
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A, 
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f 

f* 
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L 

M 
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porous  pipe  property  used  in  Equation  2-3 

pipe  cross-sectional  area  (ft*) 

2 

pipe  wall  perimeter  area  (ft  ) 
pipe  diameter  (ft) 
friction  coefficient  (Eq  3-3) 
impermeable  wall  friction  coefficient 
average  condenser  friction  coefficient  (Eq  3-16) 
axial  force  acting  on  condenser  (Eq  3-15) 


flbf~ 

L  sec 


f  lbm-  ft  ' 
conversion  factor  I - - 

(lb,-  sec*J 


pipe  length  (ft) 
Mach  number 

mass  flow  rate 


pressure 


ideal  gas  constant  for 

radial  Reynolds  number 
axial  Reynolds  number 
temperature  (°R) 
average  axial  velocity 
radial  velocity  at  wall 


(Eq  3-14) 
(Eq  3-13) 

(ft/sec) 

(ft/sec) 


v 


X 


axial  location  (ft) 


2 

A(P  )  difference  in  square  of  pressures  across  pipe  va 


AW  change  in  momentum  (lbf) 

AX  axial  increment  (ft) 

<p  momentum  flux  factor  (Eq  1-1) 

Y  ratio  of  specific  heats 

fj  dynamic  viscosity 

u  momentum  diffusivity  (ft  /sec) 

p  density 

ioVv  mass  flux 

t  wall  shear  stress 

w 

K  dimensionless  axial  velocity  (Eq  2-1) 

Superscripts  and  Subscripts ; 
ce  condenser  entrance 

d  end  of  the  condenser 

ue  end  of  the  evaporator 

u  beginning  of  the  evaporator 

t  upstream  end  of  the  AX  increment 

2  downstream  end  of  the  AX  increment 


,lbf  —  sec,. 

I  fv2  J 
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Abstract 

This  thesis  examines  the  combined  effects  of  pressure 
gradients  and  blowing  and  suction  on  frictional  forces  in  a 
heat  pipe  with  relatively  low  radial  Reynolds  numbers.  A 
porous  tube  is  used  to  simulate  the  heat  pipe  and  a  vacuum 
pump  is  used  to  generate  the  air  flow.  By  measuring  the 
static  pressure  variation  along  the  pipe  wall  and  using  a 
one-dimensional,  incompressible,  numerical  model,  the 
frictional  forces  are  obtained  and  compared  to  laminar 
fully-developed  theoretical  values.  Four  flow  rate  cases 
with  radial  Reynolds  numbers  of  1.8,  3.5,  6.5,  and  12.6  were 
studied.  In  this  range,  the  flow  in  the  evaporator  was 
fully-developed.  In  the  condenser,  however,  the 
fully-developed  solution  consistently  under  predicted  the 
average  condenser  friction  coefficient.  Deviation  from  the 
fully-developed  solution  increased  as  the  flow  rate 
increased . 
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INCOMPRESSIBLE  FLOW 
FRICTION  COEFFICIENTS 
IN  A  SIMULATED  HEAT  PIPE 

I  .  Introduction 

The  purpose  of  this  chapter  is  twofold.  First,  a 
review  of  past  efforts,  both  analytical  and  experimental, 
which  studied  flow  conditions  in  porous  constant  diameter 
circular  cylinders  with  mass  transfer  at  the  walls  will  be 
made.  This  review  will  provide  a  framework  for  under¬ 
standing  the  results  of  this  investigation  in  addition  to 
providing  sources  for  comparison.  Second,  the  objective  and 
scope  of  this  study  will  be  outlined. 

Literature  Review 

Examining  fluid  flow  in  porous  tubes  has  proven  to  be 
an  excellent  way  of  understanding  the  combined  effects  of 
pressure  gradients  and  blowing  and  suction  on  frictional 
forces  in  heat  pipes.  Most  investigators  have  concentrated 
on  the  effects  of  either  blowing  or  suction  separately.  For 
this  reason,  the  literature  discussed  will  be  separated  into 
1)  laminar  flow  in  a  porous  tube  with  blowing,  2)  laminar 
flow  in  a  porous  tube  with  suction,  and  3)  flows  simulating 
heat  pipe  operation. 
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Laminar  FI ow  in  a  Porous  Tube  With  B1 owing .  Berman  (1) 
was  the  first  to  solve  the  Navi er-St okes  equations  for 
steady,  fully-developed,  incompressible,  laminar  flow  in  a 
porous  tube  with  uniform  injection.  He  observed  that  fluid 
injection  through  the  pipe  walls  has  a  noticeable  effect  on 
veloc;.y  profiles  and  pressure  gradients.  He  showed  that 
blowing  in  all  cases  would  cause  the  pressure  to  decrease  in 
the  direction  of  axial  flow.  He  also  found  that  blowing 
would  flatten  the  velocity  profile  at  the  center  of  the  pipe 
and  increase  the  slope  of  the  profile  at  the  pipe  wall. 

Kinney  (2)  using  a  numerical  solution  to  the  boundary- 
layer  equation,  showed  that  for  fully-developed  laminar  pipe 
flow  with  either  uniform  blowing  or  suction,  the  product  of 
the  friction  coefficient,  f,  and  axial  Reynolds  number, 

Re^  is  only  a  function  of  the  radial  Reynolds  number,  Rev- 
He  termed  this  relationship  "the  universal  law  of  wall 
friction."  With  no  blowing  nor  suction,  fRe^  *  16.  He 
further  noted  that  this  product  increases  with  blowing  and 
decreases  with  suction.  Kinney  also  discussed  the 
importance  of  a  parameter  he  termed  "the  momentum  flux 
factor,"  which  is  defined  as 


This  parameter  provides  an  indication  of  the  relative 
flatness  of  the  velocity  profile.  For  the  impossible  case 


1-2 


of  slug  flow,  4>  would  equal  unity,  while  for  the  impermeable 
wall  case  in  fully-developed  laminar  pipe  flow,  4>  is  equal 
to  1.33.  For  increasing  blowing  rates,  which  have  flatter 
velocity  profiles,  <p  decreases  asymptotically  to  1.2337. 
Finally,  Kinney  verified  that  “for  a  given  axial  Reynolds 
number,  blowing  increases  the  wall  friction."  This  is 
opposite  to  that  found  for  flat-plate  boundary-layer  flows 
for  which  the  wall  friction  is  reduced  by  fluid  injection 
normal  to  the  wall.  The  reason  for  this  is  that  for 
boundary-layer  flow  over  a  flat-plate  at  zero  incidence,  the 
axial  pressure  gradient  is  zero,  even  in  the  presence  of 
uniform  blowing  normal  to  the  plate  surface.  Here,  blowing 
decreases  the  slope  of  the  velocity  profile  at  the  surface 
resulting  in  less  wall  friction.  For  the  pipe  flow  with 
blowing,  the  axial  pressure  gradient  is  negative  to  provide 
the  force  required  for  axial  acceleration.  As  the  magnitude 
of  the  pressure  gradient  increases  with  blowing,  the  slope 
of  the  velocity  profile  near  the  wall  increases.  The  end 
result  is  greater  wall  friction. 

The  first  attempt  to  verify  experimentally  the 
similarity  solutions  for  fully-developed  laminar  flow  in  a 
porous  pipe  with  wall  injection  was  performed  by  Bundy  and 
Weissberg  (3) .  They  were  able  to  confirm  the  negative  axial 
pressure  gradients  predicted  by  Berman. 

Gupta  <4)  studied  the  effect  of  blowing  on  the  laminar 
flow  in  the  inlet  region  of  a  tube.  His  interest  was  in 
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controlling  the  length  of  the  inlet  region.  Because  blowing 
increases  boundary-layer  growth,  his  proposal  was  to  use 
blowing  as  the  controlling  device.  His  analysis  assumed 
steady,  laminar  flow  of  an  incompressible  Newtonian  fluid 
entering  a  circular  tube  with  a  constant  uniform  velocity. 
The  same  fluid  was  injected  into  the  flow,  normal  to  the 
|  wall  at  a  constant  velocity.  He  found  that  as  the  blowing 

Re^  increased,  the  axial  distance  required  for  the 
attainment  of  a  fully-developed  flow  solution  decreased.  He 
also  found  that  for  a  given  blowing  Re^,  the  length  required 
to  achieve  the  fully-developed  solution  decreases  as  the 
inlet  axial  Reynolds  number  is  decreased.  This  fact  when 
applied  to  heat  pipes,  where  the  flow  begins  in  the 
evaporator  with  no  axial  velocity,  allows  the  hydrodynamic 
entrance  length  in  the  evaporator  to  be  ignored  for  even 
extremely  small  evaporation  rates. 

Bowman  (5),  determined  experimentally  that  transition 
is  retarded  due  to  mass  injection.  He  found  that  the  flow 
remains  laminar,  when  blowing  is  present,  up  to  axial 
Reynolds  numbers  on  the  order  of  10**.  Combining  Gupta’s 
observation  that,  with  blowing,  fully-developed  laminar  flow 
is  quickly  achieved  with  Bowman's  observation  that,  with 
blowing,  the  flow  remains  fully-developed  and  laminar 
provides  an  insight  as  to  why  flow  with  blowing  is 
relatively  well  behaved  compared  to  flow  with  suction. 
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Laminar  Flow  in  a  Porous  Tube  With  Suet i on .  Berman  (1) 
showed  that,  for  fully-developed  laminar  flow,  relatively 
small  amounts  of  suction  can  lead  to  a  pressure  rise  in  the 
axial  flow  direction.  In  the  presence  of  suction  at  the 
wall,  the  decrease  in  axial  momentum  tends  to  cause  an 
increase  in  pressure  in  the  flow  direction,  while  the  wall 
shear  tends  to  decrease  the  pressure.  He  demonstrated  that 
changes  in  the  familiar  parabolic  velocity  profile  observed 
in  tubes  with  impermeable  walls  are  much  more  pronounced 
when  fluid  is  removed  through  the  system  boundary  than  when 
fluid  is  injected.  For  small  suction,  the  centerline  value 
of  the  velocity  increases  and  the  slope  of  the  profile  at 
the  wall  decreases.  Because  the  velocity  of  the  fluid  near 
the  wall  is  small,  it  is  affected  more  by  the  adverse 
pressure  gradient.  As  suction  increases,  the  velocity 
gradient  at  the  wall  approaches  zero  resuming  in  a  typical 
separation  profile.  He  also  noted  that  the  amount  of 
suction  required  to  produce  a  separation  profile  is  very 
small,  (on  the  order  of  Rey  =  4) .  Moreover,  he  showed  that 
this  phenomena  is  independent  of  Rex>  In  contrast,  for 
external  flows,  where  the  pressure  gradient  is  not  dependent 
on  the  suction  rate,  suction  actually  stabilizes  the  flow. 

In  this  region  of  instability  Berman  (1)  made  his  most 
interesting  discoveries.  He  noticed  several  peculiarities 
with  his  numerical  similar  solutions.  First  when  4.1  <  Re^ 
<  4.6,  there  are  two  similar  solutions  for  the  ratio  of 
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maximum  to  average  axial  velocity,  (similar  in  the  sense 
that  the  shapes  of  the  velocity  profiles  do  not  vary  with 
distance  along  the  pipe) .  He  felt  only  one  of  the  solutions 
had  any  physical  significance  and  that  transition  to 
turbulent  flow  would  occur  before  the  second  solution  could 
ever  be  achieved.  More  importantly,  when  he  sought  similar 
solutions  for  Re^  between  the  range  of  4.6  and  about  18,  he 
found  that  none  existed.  For  Re^  greater  than  18,  he  found 
additional  multiple  solution  ranges  corresponding  to  reverse 
flow  in  various  regions  of  the  flow  between  the  wall  and  the 
centerline  of  the  pipe.  Berman  did  not  expect  the  absence 
of  solutions  for  Rew  between  4.6  and  18. 

Weissberg  (6)  was  the  first  to  investigate  the  suction 
region  identified  by  Berman  as  having  a  non-similar 
solution.  It  was  Weissberg  who  first  referred  to  this 
region  as  the  "forbidden"  range  of  suction  rates.  He  dealt 
with  the  partial  differential  equation  of  motion  for  the 
pipe  case  when  similarity  of  the  solutions  is  not  imposed. 
His  main  assumptions  were  steady,  incompressible, 
axially-symmetr ic  flow  in  a  pipe  with  uniform  suction.  He 
simplified  the  general  equation  of  motion  for  the  limit  of 
high  axial  Reynolds  number.  He  concluded  that  the 
nonexistence  of  similar  solutions  is  to  be  associated  with  a 
state  of  the  flow  field  for  which  the  entrance  region 
extends  all  the  way  from  the  pipe  inlet  to  the  location 
where  the  axial  velocity  is  reduced  to  zero  by  wall  suction. 
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Weissberg's  analysis  narrowed  this  region  to  4.8188  <  Rev  < 
15 .2688. 

Kinney  (2)  also  found  that,  for  suction  of  a  fully- 
developed  laminar  flow,  wall  friction  decreases  with 
increasing  Rew .  He  established  the  limit  for  R®w  which 
would  cause  the  wall  friction  to  be  reduced  to  zero  at  no 
less  than  4.626.  The  last  calculation  he  carried  out  was 
for  Rev  =  4.618.  This  can  be  compared  to  Berman's  original 
est imat  e  of  4.6. 

Kinney  also  provided  values  for  the  fully-developed 
laminar  <p  when  suction  is  present.  He  did  not  compute  any 
<p’ s  for  Rev  >  4.618.  For  radial  Reynolds  numbers  larger 
than  this,  fully-developed  laminar  flow  does  not  exist.  As 
suction  increases  from  zero  to  the  edge  of  the  "forbidden" 
suction  range,  he  showed  <p  for  fully-developed  laminar  flow 
increasing  from  1.33  to  approximately  1.9. 

Raithby  (7)  ,  u'sing  the  same  assumptions  as  Berman  and 
Kinney,  determined  the  range  for  which  the  wall  friction 
becomes  zero  in  fully-developed  laminar  flow  to  be  between 
Rew  of  4.5978  and  4.5980.  Raithby,  like  Berman,  had 
multiple  solutions  for  suction  radial  Reynolds  numbers 
greater  than  the  "forbidden"  suction  range.  He  felt  that 
the  velocity  profile  selected  by  the  flow  would  depend  on 
the  inlet  velocity  profile.  The  most  probable  profile  would 
likely  transition  to  turbulent  flow  due  to  a  combination  of 
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velocity  profile  inflections  and  an  opposing  pressure 
gradient . 

Hornbeck,  Rouleau,  and  Osterle  (8)  developed  a 
numerical  model  to  examine  the  entrance  region  for  a  round 
pipe  with  any  entrance  velocity  profile  and  any  suction 
boundary  condition  including  inertia  effects.  They  limited 
their  analysis  to  cases  where  the  suction  rates  were  small, 
(small  values  for  the  ratio  of  the  mean  suction  velocity  to 
the  mean  axial  velocity) .  As  in  the  previous  uniform 
suction  studies,  they  found  that  for  moderate  suction  rates, 
the  loss  of  axial  momentum  flux  resulting  from  flow  through 
the  wall  was  only  partially  offset  by  the  wall  shear  and 
the  pressure  rose  in  the  axial  direction.  For  very  small 
suction  rates,  radial  Reynolds  numbers  of  approximately  2 
and  less,  the  effect  of  wall  shear  will  be  greater  than  that 
due  to  the  loss  of  axial  momentum  flux,  and  the  pressure 
will  drop  with  axial  distance.  For  uniform  inlet 
velocities,  they  found  that  the  pressure  would  decrease 
axially  even  more  than  with  a  fully-developed  laminar 
parabolic  inlet  velocity  profile.  For  the  uniform  inlet 
velocity  profile  case,  if  Rev  is  sufficiently  large,  the 
pressure  would  first  decrease  for  some  axial  distance  and 
then  begin  to  rise.  Only  when  the  high  velocity  gradient  at 
the  wall  has  decreased  enough,  does  the  loss  of  axial 
momentum  by  flow  through  the  walls  begin  to  have  its  effect 
on  the  pressure,  causing  it  to  rise. 
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The  first  to  experiment  in  the  "forbidden"  range  of 
suction  were  Quaile  and  Levy  (9) .  They  also  performed  an 
analysis  assuming  steady,  laminar,  incompressible  flow  in  a 
constant  diameter  porous  tube  with  uniform  suction  at  the 
wall  to  compare  with  their  experimental  results.  They 
measured  the  axial  pressure  variations  and  the  velocity 
profiles  of  a  silicone  fluid  entering  a  region  of  uniform 
suction  with  a  parabolic  velocity  profile.  For  these 
conditions,  their  results  were  compatible  with  the  earlier 
theories.  Provided  in  their  report  is  a  figure  showing 
theoretical  axial  locations  where  separation  should  occur 
for  a  given  Rew .  For  Rev  between  6  and  22,  the  measured 
separation  locations  occurred  slightly  later  than  predicted. 
For  Re^  =  30,  the  measured  separation  location  occurred 
slightly  earlier  than  predicted. 

Flows  Simulating  Heat  Pipe  Qperat i on .  Measurement  of 
vapor  dynamics  in  actual  heat  pipes  presents  many 
experimental  difficulties.  Because  heat  pipes  are  carefully 
sealed  from  their  environment,  it  is  difficult  to  determine 
the  velocity  profiles  in  the  vapor  region.  Wageman  and 
Guevara  (10)  were  able  to  demonstrate  that  injecting  a  fluid 
through  the  walls  over  half  the  length  of  a  porous  pipe  and 
then  extracting  the  fluid  through  the  walls  over  the  second 
half  of  the  porous  pipe  simulated  the  vapor  dynamics  in  a 
heat  pipe  quite  well.  The  mass  injection  at  the  wall 
simulates  the  evaporator  and  the  mass  extraction  at  the  wall 
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simulates  the  condenser.  For  heat  pipes,  the  velocity 
profile  leaving  the  evaporator  is  somewhere  between  the 
parabolic  and  uniform  velocity  profiles  and  unless  the 
adiabatic  region  is  sufficiently  long,  the  parabolic  profile 
is  never  attained  at  the  inlet  of  the  condenser.  Since 
Wageman  and  Guevara,  other  researchers  such  as  Quaile  and 
Levy,  and  Bowman  have  used  a  porous  pipe  to  simulate  a  heat 
pipe  . 

Kinney  (2)  as  mentioned  previously,  developed  a 
relationship  for  the  friction  coefficient  and  axial  Reynolds 
number  product  for  various  radial  Reynolds  numbers  for 
fully-developed  laminar  flow.  Bowman  (5)  modified  this 
incompressible  relationship,  based  on  experimental  data,  to 
extend  to  compressible  flow.  For  laminar  flow,  he  found  the 
relationship  to  be 
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For  turbulent  flow  with  axial  Reynolds  numbers  less  than 
150,000  Kinney  and  Sparrow  (11)  suggested  the  following 
relationship 
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where  V  is  the  radial  velocity  at  the  wall,  U  is  the 
w 

* 

average  axial  velocity,  and  f  ,  the  impermeable  wall 
friction  coefficient,  may  be  estimated  by  the  Blasius 
equation 

0.079 

f*  =  -  (1-4) 


The  primary  objective  of  this  research  is  to  study  the 
behavior  of  a  simulated  heat  pipe  operating  at  relatively 
low  radial  Reynolds  numbers.  To  better  understand  this 
behavior,  data  from  a  simulated  heat  pipe  will  be  processed 
with  a  one-dimensional,  numerical  model.  To  lower  the 
radial  Reynolds  number  via  a  lower  density,  air  will  be 
pulled  through  the  porous  tube  with  a  vacuum  pump  instead  of 
using  compressed  air  upstream.  The  static  pressure  will  be 
measured  along  the  pipe  wall.  This  pressure  distribution 
becomes  the  input  to  a  one-dimensional,  incompressible, 
numerical  model  for  a  constant  diameter  porous  tube  with 
variable  mass  injection  and  extraction.  The  numerical  model 
will  provide  the  friction  coefficient  distribution  in  the 
evaporator  as  well  as  an  average  friction  factor  for  the 
condenser.  These  friction  coefficients  will  then  be 
compared  to  the  fully-developed  solutions  predicted  by 
Kinney.  Also,  the  measured  axial  separation  locations  will 
be  compared  to  those  predicted  by  Quaile  and  Levy. 
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II .  Experimental  Model 


The  purpose  of  this  study  was  to  obtain  static  pressure 
data  (axial  pressure  distribution)  for  flow  inside  a  porous 
tube  simulating  a  heat  pipe  at  relatively  low  radial 
Reynolds  numbers.  The  pressure  measurements  were  used  as 
input  to  a  one-dimensional  numerical  model  discussed  in 
detail  in  Chapter  3.  In  this  chapter,  a  description  of  the 
porous  pipe,  a  discussion  of  the  experimental  set  up,  and 
the  calibration  results  will  be  presented. 

Por ous  Pipe 

Following  the  example  of  Bowman  (5) ,  an  ultra  high 
molecular  weight  polyethylene  porous  pipe  was  obtained  from 
Porex  Technologies ,  (PVC  tube  part  #  5139).  The  pipe  had  an 
inside  diameter  of  0.5  inches  and  a  wall  thickness  of  0,25 
inches.  The  average  pore  size  was  20  microns.  The  length 
used  in  the  heat  pipe  simulation  was  31.81  inches.  The 
evaporator  and  condenser  sections  were  both  15.66  inches 
long  separated  by  an  adiabatic,  impermeable  wall  of  0.50 
inches  in  length.  Both  ends  were  sealed  with  plastic  caps 
to  insure  that  all  flow  entering  or  leaving  the  pipe  was 
through  the  pipe  walls.  Half  of  the  pipe  was  then  placed 
into  a  vacuum  tank,  with  the  opposite  hall  exposed  to 
ambient  air  (see  Figure  2-1) .  A  portable  shop  vacuum  pump 
connected  to  the  vacuum  tank  generated  the  flow. 
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Pressure  Taps 

Thirty  pressure  taps  were  installed  in  the  porous  pipe 
by  cementing  0.5  inch  long  0.060  inch  outside  diameter 
stainless  steel  tubes  into  0.058  inch  outside  diameter  holes 
drilled  in  the  tube.  This  produced  a  snug  fit.  The 
stainless  steel  tube  was  inserted  in  the  porous  pipe  so  that 
its  end  was  flush  with  the  inside  wall  of  the  pipe.  Two  of 
the  thirty  pressure  taps  were  installed  through  the  end  caps 
to  obtain  pressure  measurements  for  Zf  =  0.0  and  Zf  =  1,0 
where  the  dimensionless  axial  location,  is  defined  as 


The  remaining  twenty  eight  pressure  taps  were  distributed 
axially  along  the  pipe,  eleven  in  the  evaporator  and 
seventeen  in  the  condenser  (see  Figure  2-2) ,  with  their 
axial  locations  documented  in  the  numerical  code  found  in 
Appendix  B. 

Pressure  Measurement s 

All  pressure  measurements  were  made  using  a  T-type 
Scanivalve,  from  Scanivalve  Corporation.  A  Robinson  Halpern 
model  157B-WQ20D-F-V31  differential  pressure  transducer  with 
a  pressure  range  of  +/-  0-2  inches  of  water  producing  an 
output  of  +/-  1  volt  was  used  in  the  Scanivalve.  The  system 
allowed  for  36  pressure  measurements  with  the  one  pressure 
transducer.  Two  of  the  six  remaining  scanivalve  ports  were 
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used  to  measure  the  pressure  in  the  vacuum  tank,  and  the 
ambient  pressure.  For  the  three  lowest  flow  rate  cases,  all 
pressure  differences  were  measured  with  respect  to  ambient 
pressure.  For  the  highest  flow  rate  case,  all  pressure 
differences  were  measured  with  respect  to  the  pressure  port 
immediately  leaving  the  evaporator.  To  obtain  absolute 
pressures,  the  ambient  pressure  was  measured  with  a 
barometer.  The  main  disadvantage  of  the  system  was  the 
amount  of  time  required  by  the  transducer's  diaphragm  to 
reach  a  steady  pressure  reading.  The  pressure  transducer 
was  calibrated  using  an  inclined  water  manometer  as  the 
standard . 

Data  Acguisit i on 

Data  acquisition  was  done  via  a  Zenith  Z-100  computer. 
Pressure  transducer  input  into  the  computer  was  digitized 
through  a  Dual  Systems  Control  model  AIM12  analog-to-digital 
card  in  the  computer.  An  offset  option  on  the  transducer 
and  the  card's  amplification  were  set  so  signals  from  0  to 
+/-  1  volt 8  could  be  read  with  a  resolution  of  0.04883 
millivolts.  Shielded  cables  were  used  to  reduce  noise 
introduced  into  the  system.  Noise  was  further  reduced  by 
averaging  1000  data  samples  for  every  pressure  reading.  To 
allow  time  for  the  transducer's  diaphragm  to  reach 
equilibrium,  a  delay  loop  was  used  as  the  scanivalve  changed 
ports.  Once  the  data  had  been  read  and  averaged  by  the 
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computer,  it  was  displayed  on  the  monitor  and  stored  on  an 
electronic  disk. 

Pipe  Cal ibrat i on 

In  addition  to  the  pressure  distribution  along  the 
porous  pipe,  the  numerical  model  required  information 
regarding  the  resistance  of  the  porous  pipe  to  mass  flux 
(pV  )  through  its  walls.  Muskat  (12)  showed  using  Darcy’s 
Law  that,  for  the  flow  of  a  compressible  gas  through  a 
porous  medium,  pV^  can  be  related  to  the  pressure  difference 
across  the  porous  medium  by  the  expression 

A(P2)  =  A(pVw)2  +  B(pVv)  (2-2) 

where  the  constants  A  and  B  are  properties  of  the  medium. 

For  all  cases  of  this  study,  the  mass  flux  rates  were  very 
small  and  the  contribution  of  the  A(pV^)  term  was 
negligible.  Muskat’s  relationship  was  therefore  simplified 
t  o 

A(P2)  =  B(pVw)  (2-3) 

To  determine  the  constant  B,  A(P2)  and  pVw  were 
measured  at  various  mass  flow  rates  .  The  two  2  inch  samples 
used  during  the  calibration  were  obtained  from  opposite  ends 
of  the  porous  pipe  used  to  simulate  the  heat  pipe.  Two 
pressure  ports  for  each  sample  connected  to  a  50  inch  U-tube 
water  manometer  were  used  to  measure  the  average  pressure 
difference  across  the  sample  walls.  A  1  inch  precision  bore 
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"flovrator"  with  a  range  of  0.05  -  0.44  cfm  was  used  to 

2 

measure  pV  .  The  method  used  to  measure  £(P  )  and  pV  was 

V  v 

the  same  as  outlined  by  Bowman  (5)  with  the  exception  that 
flow  was  induced  by  a  vacuum  pump  instead  of  compressed  air. 
It  was  assumed  that  the  temperature  of  the  air  entering  the 
"flowrator"  was  the  same  as  ambient.  Using  a  least-squares 
technique  to  curve  fit  the  experimental  data  shown  in  Figure 
2-3,  the  constant  B  in  Equation  (2-2)  was  found  to  be 

3.592472958  X  10°  - •  As  can  be  seen  in  Figure  2-3,  the 

bob  indicator  tended  to  stick  to  the  walls  of  the 
"flowrator"  at  the  lower  flow  rates.  When  calibrating 
sample  1,  readings  were  taken  from  the  lowest  to  highest 
flow  rates.  When  calibrating  sample  2,  readings  were 
reversed . 
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Figure  2—3.  Porous  Pipe  Calibration  Results. 


III.  Numerical  Model 


This  chapter  describes  the  one-dimensional  numerical 
model  used  to  reduce  the  data  measured  during  the 
experimental  phase  of  the  study.  First,  the  governing 
equations  used  in  solving  the  flow  properties  will  be 
presented.  Then,  the  solution  technique  will  be  described. 
Governing  Eauat ions 

One-dimensional,  steady,  adiabatic,  incompressible  flow 
was  assumed  for  the  model.  In  the  numerical  simulation  code 
contained  in  Appendix  B,  a  provision  for  compressible  flow 
has  been  made  following  the  example  of  Holladay  (13) .  All 
flows  evaluated  in  this  study,  however,  are  incompressible. 

To  approximate  the  shear  stress  at  the  pipe  wall, 
Newton’s  second  law  was  applied  to  an  element  of  fluid  in  the 
pipe  as  shown  in  Figure  3-1 . 


Figure  3-1.  Element  of  Fluid  for  Equation  of  Motion 

Applying  Newton’s  second  law  for  steady  flow  to  the  fluid 
element  gives 
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where  A_  is  the  wall  surface  area  for  the  increment  AX, 
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Solving  for  gives 


T  = 

V 


m  U  <p 

1  11 


"«UA 


(p.  '  PJ‘ 


(3-2) 


where  A  is  the  cross-sectional  area  equal  to  A  and  A  .  The 

C  12 

friction  coefficient  is  given  by 
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where 


M 


(3-4) 


Substituting  the  expression  for  given  by  Equation  (3-2)  in 
Equation  (3-3)  gives  an  expression  for  the  friction  coefficient 
in  terms  of  quantities  which  can  be  calculated  from  the  input 
pressure  distribution 
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Because  4>  in  the  evaporator  was  known  and  nearly 
constant,  the  total  force  for  the  condenser  could  be 


det  ermined . 
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where  the  subscripts  u,  »,  and  ce  represent  the  beginning  of 
the  evaporator,  the  end  of  the  condenser,  and  the  condenser 
entrance,  respectively.  By  using  average  condenser  values 
for  axial  velocity  and  density,  the  average  friction 
coefficient  for  the  condenser  could  be  evaluated  using 
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The  average  condenser  axial  velocity  was  calculated  as  half 
the  axial  velocity  entering  the  condenser. 


Sol ut i on  Method 

A  marching  technique  was  used  to  solve  for  flow 
properties  along  the  simulated  heat  pipe.  The  static 
pressure  data  obtained  from  the  experimental  work  was  first 
read  into  the  computer  code.  Next,  the  program  prompted  for 
average  values  of  <p  for  the  blowing  and  suction  regions. 

With  predefined  constants  for  the  pipe  geometry  and  air 
properties  available,  and  knowing  the  initial  conditions  at 
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£  *  0.0,  the  set  of  governing  equations  were  used  to  solve 
for  the  downstream  Mach  number,  axial  Reynolds  number, 
radial  Reynolds  number,  and  friction  coefficient  for  each  AX 
increment  along  the  pipe. 

Three  subroutines  were  utilized  to  1)  interpolate  for 
intermediate  pressures,  and  solve  for  flow  properties  for 
either  2)  incompressible,  or  3)  compressible  flow.  These 
subroutines  were  individually  verified  using  known  input  and 
output  values.  The  main  program  was  verified  by  using  the 
pure  mass  case  (no  friction).  Using  Shapiro's  influence 
coefficients  for  pure  mass  addition  (14),  with  constant 
specific  heat  and  molecular  weight 


M-/2[r+1) 

m  =  m  -  (3-17) 

l+j'M 

and 

*  r+1 

P  =  P  -  (3-18) 

1+rM2 


and  for  a  constant  wall  mass  flux,  m 


dm  _  dx 

•  x 

in 

the  Mach  number  as  a  function  of  axial  location, 
found  using 
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x,  can  be 
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Selecting  the  Mach  number  at  the  end  of  the  evaporator,  le , 

* 

determines  x  and  specifies  a  Mach  number  for  each  axial 


location.  The  pressure  variation  is  then  known  using 


P 

P  =  -  (3-21) 

1+yM 


where  Po  is  the  pressure  at  the  end  of  the  tube.  Solving  for 

.  I  • 

the  constant  mass  flux,  m  ,  gives 


(3-22) 


(3-23) 


The  static  pressure  in  the  condenser  for  the  pure  mass  case 
mirrors  the  evaporator  solution. 

For  the  experimental  cases,  the  numerical  code  provided 
two  additional  checks.  First,  the  force  balance  over  the 
entire  pipe  was  calculated.  When  considering  the  pipe 
system,  the  only  forces  acting  on  the  fluid  are  the 
pressures  at  the  ends  of  the  pipe  and  the  total  wall  shear 
force . 
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D 

Force  Balance  =  [pw  -  Pd]ac  -  J  (TwAe).  (3'24) 

1=0 

Ideally,  the  Force  Balance  should  be  zero  since  the  drop  in 
pressure  between  the  two  ends  of  the  pipe  is  a  result  of  the 
wall  shear  stress.  A  relative  percent  error  was  provided  by 
dividing  the  Force  Balance  by  the  pressure  force  and 
multiplying  by  100.  The  final  check  made  was  a  mass  balance 
over  the  entire  pipe  to  see  if  the  total  mass  injected  into 
the  pipe  equaled  the  total  mass  removed  from  it.  A  relative 
percent  error  was  provided  by  dividing  the  mass  balance  by 
the  mass  entering  the  condenser  and  multiplying  by  100. 
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IV .  Discussion  of  Results 

As  can  be  seen  from  Figure  4-1,  Rey  for  all  four  test 
cases  was  relatively  constant.  As  the  flow  rates  increased, 
the  axial  variation  of  Rev  became  more  obvious.  The  largest 
flow  rate  case  had  a  0.9  %  change  in  Rew  between  the 
beginning  of  the  evaporator  and  the  end  of  the  condenser. 

Figure  4-2  shows  how  the  pressure  varied  axially  along 
the  walls  of  the  simulated  heat  pipe  for  all  four  Rew  cases. 
The  data  used  to  create  Figure  4-2  can  be  found  in  Appendix 
A.  The  pressure  variation  obtained  along  the  pipe  was  a 
result  of  two  effects:  1)  the  friction  encountered  along 
the  pipe  walls,  and  2)  the  rate  of  change  in  momentum  along 
the  pipe.  For  flow  without  friction,  the  pressure  would 
drop  as  the  flow  accelerated  in  the  evaporator  and  would 
completely  recover  as  the  flow  decelerated  in  the  condenser. 
The  total  effect  of  the  wall  friction  force  can  be  seen  as  a 
pressure  deficit  at  the  end  of  the  pipe.  In  each  of  the 
four  Rew  cases,  the  flow  did  not  achieve  complete  pressure 
recovery.  The  pressure  drop  due  to  friction  was  found  to 
increase  with  increased  mass  flow  rate.  For  the  lowest  flow 
rate  case,  (Rew  =  1.8),  the  pressure  in  the  condenser 
dropped  instead  of  rising.  This  occurred,  as  predicted  by 
Hornbeck,  et .  al .  for  suction  radial  Reynolds  numbers  less 
than  approximately  2  (8) ,  because  the  effect  of  wall  shear 
is  greater  than  that  due  to  the  loss  of  axial  momentum  flux. 
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Figure  4—2.  Experimental  variation  of  wall  static 
pressure  with  axiai  location  along 
porous  pipe. 


As  can  be  seen  in  Figure  4-2,  flow  separation  occurred 
in  the  condenser  for  the  Rev  of  6.5  and  12.6  cases.  To 
compare  separation  locations  with  those  predicted  and 
measured  by  Quaile  and  Levy  (9),  their  non-dimensional 
coordinate  system  had  to  be  converted  to  that  used  in  this 
experiment  since  they  did  not  have  an  evaporator  region 
preceding  their  suction  region.  Figure  4-3  has  three  pieces 
of  information:  1)  a  curve  fit  of  Quaile  and  Levy's 
converted  theoretical  separation  predictions,  2)  the 
converted  experimental  separation  locations  found  by  Quaile 
and  Levy,  and  3)  the  separation  locations  found  in  this 
experiment.  Recall  that  Quaile  and  Levy  assumed  uniform 
suction  and  a  parabolic  inlet  velocity  in  their  derivation 
of  theoretical  separation  locations.  The  suction  in  this 
experiment  was  nearly  uniform  as  can  be  seen  in  Figure  4-1, 
however,  since  the  flow  entering  the  suction  region  was 
preceded  by  a  blowing  region  and  a  short  impermeable  wall 
with  an  L/D  of  1,  the  velocity  profile  entering  the 
condenser  will  be  flatter  than  parabolic,  <  1.33)  .  The 
velocity  profile  entering  the  suction  region  requires  some 
axial  distance  to  adjust  to  local  wall  conditions.  This 
would  explain  why  separation  in  this  experiment  occurred 
slightly  after  that  predicted  by  Quaile  and  Levy.  The  two 
lowest  flow  rate  cases  had  Re^’s  less  than  4.6,  the 
approximate  Rev  which  would  cause  separation  in  a 
fully-developed  flow.  No  flow  separation  occurred  for  these 
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Figure  4—3.  Curve  Ft  of  Quaile  and  Levy's  Prediction 
of  Separation  for  constant  wail  suction 
and  parabolic  inlet  velocity  profile. 
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Friction  Coefficient  From  Equation  (1—2) 


1— D  Friction  Coefficient  From  Equation  (3—5) 


Figure  4—4.  Evaporator  (Blowing)  Friction  Coefficient 
Results. 


cases.  For  the  two  cases  where  separation  occurred,  the 
axial  location  plotted  in  Figure  4-3  corresponds  to  the 
center  of  the  two  pressure  ports  detecting  separation. 

Figure  4-4  is  a  graph  of  the  evaporator  friction 
coefficients  calculated  comparing  the  expression  for 
fully-developed  laminar  flow.  Equation  (1-2) ,  with  the 
one-dimensional  expression.  Equation  (3-5).  For  all  four 
cases,  the  flow  was  fully-developed  within  four  AX 
increments  along  the  pipe.  It  should  be  noted,  that  the 
values  calculated  using  Equation  (3-5)  are  highly  dependent 
on  <f> .  The  values  of  <p  used  were  obtained  from  a  graph 
Kinney  (2)  produced  numerically  assuming  fully-developed 
laminar  flow  with  constant  blowing  Re^’s.  Bowman’s  Equation 
(1-2)  is  an  empirical  curve-fit  of  Kinney's  results.  The 
values  obtained  from  Equation  (3-5)  differ  slightly  from 
those  predicted  by  Kinney.  However,  by  choosing  a  more 
appropriate  4> ,  Kinney’s  results  can  be  obtained.  The  first 
data  points  for  each  Rey  case  plotted  correspond  to  the 
third  increment  of  the  numerical  data  reduction  program. 
These  Tour  data  points  deviate  the  farthest  from  the  line 
for  fully-developed  laminar  flow  when  compared  to  respective 
Rev  data  points 

Figure  4-5  is  a  graph  comparing  Kinney's  (2) 
theoretical  universal  law  of  wall  friction  for 
fully-developed  laminar  flow  in  porous  tubes  with  uniform 
suction,  with  the  average  condenser  friction  coefficients, 
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Figure  4—5  Curve  fit  of  Kinney's  universal  law  of  wall 

friction  for  fully— developed  laminar  flow  in 
porous  tubes  and  average  condenser  values 
from  this  experiment. 


f  ,  calculated  using  Equation  (3-16) .  For  the  two  lowest 
flow  rate  cases,  f  is  slightly  higher  than  that  predicted 
by  Kinney.  This  can  be  explained  by  the  fact  that  the  flow 
entering  the  condenser,  due  to  the  blowing  in  the 
evaporator,  has  a  flatter  velocity  profile  than  the 
parabolic  inlet  velocity  profile  assumed  by  Kinney.  During 
the  length  of  the  condenser  where  the  flow  is  adjusting  to  a 
fully-developed  flow,  the  friction  should  be  higher  than  the 
fully-developed  values.  For  the  last  two  flow  rate  cases, 
higher  friction  due  to  the  adjusting  velocity  profiles,  in 
addition  to  separated  flow  would  account  for  the  larger 
deviations  from  Kinney's  theoretical  curve.  The  large 
friction  coefficients  after  separation  would  more  than 
balance  the  small  friction  coefficients  in  the  regions 
preceding  separation. 

Appendix  C  contains  printouts  for  the  four  flow  rates. 
These  printouts  indicate  the  average  values  of  <p  used  in  the 
blowing  and  suction  regions.  They  also  show  the  relative 
percent  errors  for  the  force  and  mass  balances  encountered 
during  the  experimental  and  numerical  data  reduction  phases. 
The  maximum  error  was  0.0013  %  in  the  force  balance  and 
0.0312  %  in  the  mass  balance. 
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V .  Conclusions  and  Recommendations 


Conclusions 

The  results  obtained  from  this  investigation  provide 
insights  into  the  operation  of  heat  pipes  at  relatively  low 
radial  Reynolds  numbers .  The  fully-developed  laminar  flow 
solution  can  be  used  in  the  evaporator  to  obtain  friction 
coef f icients .  In  the  condenser,  however,  the  fully- 
developed  solutions  may  not  accurately  predict  wall 
friction.  Above  radial  Reynolds  numbers  of  4.6,  the 
hydrodynamic  entrance  regions  can  be  a  major  portion  of  the 
condenser  length.  In  heat  pipes,  unless  there  exists  a 
sufficiently  long  adiabatic  region,  the  flow  entering  the 
condenser  has  flatter  velocity  profiles  than  those  in  the 
impermeable  wall  case.  Through  the  length  of  the  condenser 
where  the  flow  is  adjusting  to  local  wall  conditions,  the 
wall  friction  will  be  higher  than  the  fully-developed 
solution.  Also,  because  of  this  phenomenon,  flow  separation 
will  be  slightly  delayed.  In  this  investigation,  an  average 
condenser  friction  factor  was  found  for  radial  Reynolds 
numbers  less  than  13.  This  result  should  prove  helpful  for 
design  purposes. 

Recommendations 

Some  ideas  for  further  research  are: 

1)  Experimentally  determine  condenser  flow  separation 
locations  and  average  friction  coefficients  for  radial 
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Reynolds  numbers  greater  than  13,  where  blowing  and  suction 
rates  are  not  constant,  and  the  condenser  inlet  velocity 
profile  becomes  flatter. 

2)  Determine  the  onset  of  transition  to  turbulent  flow 
and  its  effect  on  average  condenser  friction  coefficients. 

3)  Perform  a  numerical,  two-dimensional  flow  analysis, 
using  experimental  pressure  data,  to  determine  appropriate 
momentum  flux  factors,  <p ,  for  the  developing  flow  in  the 
condenser . 


i 
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Appendix  A:  Experimental  Data 
This  appendix  contains  raw  data  for  the  four  flow  rate 
cases  obtained  during  the  experimental  portion  of  this 
study.  Test  run  numbers  1,  2,  3,  and  4  correspond  to  radial 
Reynolds  numbers  of  3.5,  1.8,  6.5,  and  12.6  respectively. 

All  pressures  are  given  in  psia  and  all  temperatures  are 
given  in  °R .  The  following  format  is  used  for  all  four  test 


cases  : 


test  run  number 
atmospheric  pressure 
ambient  temperature 
vacuum  tank  pressure 


port 

1  1 

pressure,  ? 

=  0.0 

port 

2  pressure 

port 

3  pressure 

port 

4  pressure 

port 

5  pressure 

port 

6  pressure 

port 

7  pressure 

port 

8  pressure 

port 

9  ] 

pressure 

port 

10 

pressure 

port 

11 

pressure 

port 

12 

pressure , 

(last 

port 

13 

pressure , 

(first 

port 

14 

pressure 

port 

15 

pressure 

port 

16 

pressure 

port 

17 

pressure 

port 

18 

pressure 

port 

19 

pressure 

port 

20 

pressure 

port 

21 

pressure 

port 

22 

pressure 

port 

23 

pressure 

port 

24 

pressure 

port 

25 

pressure 

port 

26 

pressure 

port 

27 

pressure 

port 

28 

pressure 

port 

29 

pressure 

port 

30 

pressure , 

*  *  1. 

(beginning  of  evaporator) 


port  in  evaporator) 
port  in  condenser) 


(end  of  condenser) 


A-l 


1 

14.282667000000 

526.50 

14 .244351480722 
14 .263613999968 
14 .263609987649 
14.263602540310 
14 .263591359022 
14 .263576271017 
14 .263557421172 
14  .263534520094 
14  .263528299893 
14 .263521720020 
14.263514807101 
14.263507097221 
14  .263497110700 
14.263497081109 
14.263498061892 
14.263499050988 
14 .263500745013 
14  .263504779938 
14  .263507824992 
14  .263509883195 
14  .263511303314 
14  .263512510296 
14  .263513600389 
14.263514573218 
14.263515434911 
14.263516359503 
14  .263517049930 
14.263517850391 
14  .263518084766 
14.263518090997 
14  .263518091370 
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14  .248416000000 
526.00 

14  228081934722 
14  238297320471 

14  .238295723592 

14 .238292754917 
14  .238288259536 
14  238282241124 
14 .238274679967 
14  238265590919 
14  .238263080953 
14  238260470171 
14.238257765012 
14  .238254963998 
14  238251650102 
14  238250005044 
14  .238249248215 
14 ,238248648028 
iV  238248065997 
14  .238246948099 
14  238245915972 
14  .238244968116 
14  238244105078 
14  .238243338026 
14  238242661110 
14  238242065994 
14  238241548503 
14  .238241108511 
14 .238240752108 
14  .238240482496 
14 .238240294616 
14  .238240190500 
14 .238240163414 


3 

14.302239000000 

526.00 

14 .230802744722 
14 .266756940340 
14 .266746082953 
14.266726010212 
14 .266695605500 
14 .266654991066 
14  .266604022328 
14.266542640400 
14 .266525647021 
14.266508007998 
14.266489624021 
14 .266470555221 
14  .266444511202 
14 .266431792928 
14 .266445921015 
14 .266455020220 
14.266465723482 
14 .266485780049 
14  .266503172131 
14 .266521909861 
14  .266537951121 
14.266552664012 
14  .266565232023 
14.266576731019 
14 .266586102102 
14.266597460074 
14  .266605711573 
14 .266611421120 
14  .266603292385 
14 .266604100374 
14.266604363487 


r 


4 

14.272881000000 

526.00 

14.133008654431 
14.203597320072 
14.203562960836 
14 .203499542432 
14  .203403802361 
14.203275260128 
14 .203113651026 
14  .202918830572 
14  .202864878298 
14  .202808843901 
14  .202750764670 
14.202690523947 
14 .202602026204 
14  .202584598832 
14.202647510237 
14  .202695511205 
14  .202740492984 
14.202826570110 
14  .202908462280 
14  .202985260293 
14.203056868943 
14  .203121392452 
14  .203178934831 
14  .203229592391 
14.203273830120 
14  .203274341038 
14  .203291310250 
14  .203305581193 
14 .203318130044 
14  .203325605621 
14 .203327728798 
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D  Line# 

1 

2  C 

3  C 

4  C 

5  C 

6  C 

7  C 

8  C 

9  C 

10  C 

11  C 

12  C 

13  C 

14  C 

15  C 

16  C 

17  C 

18  C 

19  C 

20  C 

21  C 

22  C 

23  C 

24  C 

25  C 

26  C 

27  C 

28  C 

29  C 

30  C 

31  C 

32  C 

33  C 


>endix  B 


10-19-88 

1  7  Microsoft  F0RTRAN77  V3.20  02/84 

PROGRAM  ONED 

TITLE:  ONE-DIMENSIONAL  HEAT  PIPE  SIMULATION 


BY :  Captain  David  Manley 
MS  Thesis,  1988 

This  program  uses  one-dimensional  incompressible 
and  compressible 

flow  models  to  solve  for  the  flow  properties  and 
friction 

coefficients  in  a  closed  porous  pipe  with  blowing 
and  suction. 

The  incompressible  model  is  used  to  calculate  the 
flow  properties 

until  the  Mach  number  reaches  0.01.  Choosing  0.01 
was  a  purely 

arbitrary  decision.  For  Mach  numbers  above  0.01, 
the  compressible 

model  is  used.  The  compressible  model  utilizes 
Shapiro’ s  method 

of  influence  coefficients  to  solve  for  the  flow 
properties . 


*****  DEFINITION  OF  VARIABLES  ***** 


A 


AC 

AP 

BETA 

CL 

CC1 


C2 

D 

DMOM 

DPS 


Constant  obtained  from  pipe 
calibration  used  in 
determining  wall  mass  flux 
(lbf/f t/sec) 

Cross-sectional  area  of  the  pipe 
(ft2) 

Surface  area  of  DX  increment  (ft2) 
Ratio  of  radial  Reynolds  number  to 
axial  Reynolds  number 
Condenser  Length  (ft) 

Intermediate  value  used  during 
pressu&  e 

interpolation  subroutine  (CC2  and 

CC3  are 

similar) 

Speed  of  sound  squared  (ft2/sec2) 
Pipe  inside  diameter  (ft) 

Change  in  momentum  across  DX 
increment  (lbf) 

(Pambient) **2  -  (P2)**2 
(lbf /ft2) **2 
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34 

C 

DX 

Spacial  step  sise  (ft) 

35 

C 

DW 

Change  in  mass  flow  rate  due  to 
mass  transfer  at 

36 

c 

the  wall  (lbm/sec) 

37 

c 

DWSUM 

Summation  of  all  the  incremental 
mass  flow  rate 

38 

c 

changes  over  the  entire  pipe  due 
to  mass  transfer 

39 

c 

at  the  wall  (lbm/sec) 

40 

c 

F 

Friction  factor 

41 

c 

FBARC 

Average  condenser  friction  factor 

42 

c 

FBL 

Bowman's  friction  factor  for 
fully-developed 

43 

c 

laminar  flow 

44 

c 

FBLREX 

Product  of  FBL  and  axial  Reynolds 
number 

45 

c 

FBT 

Kinney’s  friction  factor  for 
fully -developed 

46 

c 

turbulent  flow 

47 

c 

FBTREX 

Product  of  FBT  and  axial  Reynolds 
number 

48 

c 

FCREXC 

Product  of  FBARC  and  REXBRC 

49 

c 

FSTAR 

No  blowing/suction  fully-developed 
turbulent 

50 

c 

friction  coefficient,  (Blasius 
equation) 

51 

c 

FUD 

Constant  to  account  for  the 
statistical  variation 

52 

c 

of  the  porous  pipe  calibration 
constant  A 

53 

c 

GAM 

Ratio  of  specific  heats 

54 

c 

GC 

Newton's  gravitational  constant 
(lbm-f t/lbf /sec2) 

55 

c 

Ml 

Mach  number  at  upstream  end  of  DX 
increment 

56 

c 

M2 

Mach  number  at  downstream  end  of 

DX  increment 

57 

c 

MIS 

Ml  squared 

58 

c 

M2S 

M2  squared 

59 

c 

MFF 

Momentum  flux  factor,  (see  Kinney 
pg  1398 ; I JHMT 

60 

c 

Vol  11  1968)  ,  An  indication  of 
the  velocity 

61 

c 

profile  or  shape 

62 

c 

MFF1 

Momentum  flux  factor  at  upstream 
end  of  DX 

63 

c 

increment 

64 

c 

MFF2 

Momentum  flux  factor  at  downstream 
end  of  DX 

65 

c 

increment 

66 

c 

MFFB 

Blowing  momentum  flux  factor 

67 

c 

MFFS 

Suction  momentum  flux  factor 
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68 

C 

MS  BAR 

Average  of  MIS  and  M2S 

69 

C 

N 

Test  Run  number 

70 

C 

PATM 

Atmospheric  pressure  Centered  in 
psia  then 

71 

C 

converted  to  lbf/ft2) 

72 

C 

PEF 

Percent  error  for  force  balance 

73 

C 

PEM 

Percent  error  for  mass  balance 

74 

c 

PF 

Pressure  force  on  entire  pipe  due 
to  the 

75 

c 

difference  of  pipe  end  pressures 

76 

c 

PL 

Pipe  length  (ft) 

77 

c 

PND (1) -PND (30) 

Non-dimensional  static  pressure, 
P2/P  (X  =  0 , 0) 

78 

c 

PI 

Static  pressure  at  upstream  end  of 
DX  increment 

79 

c 

Clbf/ft2) 

80 

c 

P2 

Static  pressure  at  downstream  end 
of  DX  increment 

81 

c 

(lbf/f t2) 

82 

c 

P (1)  -P (30) 

Pressure  measurements  at 
corresponding  axial 

83 

c 

positions,  X(l)-X(29)  (entered  in 
psia  and  then 

84 

c 

converted  to  lbf/ft2) 

85 

c 

R 

Ideal  gas  constant  for  air 
(lbf-ft/lbm/R) 

86 

c 

REW 

Radial  Reynolds  number  based  on 
pipe  inside 

87 

c 

diameter 

88 

c 

REX 

Axial  Reynolds  number  based  on 
pipe  inside 

89 

c 

diamet  er 

90 

c 

REXBRC 

Average  axial  Reynolds  number  for 
the  condenser 

91 

c 

based  on  pipe  inside  diameter. 

REX  entering 

92 

c 

condenser  divided  by  2 

93 

c 

RHOC 

Representative  density  of  air  for 
condenser 

94 

c 

(lbm/ft3) 

95 

c 

RHOl 

Density  of  air  at  beginning  of  DX 
increment 

96 

c 

(lbm/ft3) 

97 

c 

RH02 

Density  of  air  at  end  of  DX 
increment  (lbm/ft3) 

98 

c 

RHOV1 

Mass  flux  through  the  wall  at  the 
upstream  end  of 

99 

c 

the  DX  increment  (lbm/f t2/sec) 

100  c 

RHOV2 

Mass  flux  through  the  wall  at  the 
downstream  end 

101 

.  c 

of  the  DX  increment  (lbm/ft2/sec) 

102  C 

RHOVAV 

Average  mass  flux  for  DX 
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103  C 


RHOBAR 


104  C 

105  C 

106  C 

107  C 

108  C 

109  C 

110  C 

111  C 

112  C 

113  C 

114  C 

115  C 

116  C 

117  C 

118  C 

119  C 

120  C 

121  C 


122  C 

123  C 

124  C 

125  C 

126  C 

127  C 

128  C 

129  C 

130  C 

131  C 

132  C 

133  C 

134  C 

135  C 

136  C 

137 


RMU 

SF 

SFC 

SFEA 

SFSUM 

T 

TEST 

TW 

TWBARC 

TO 

UBARC 

U1 

U2 

W1 

W2 

WBAR 

WCK 

XL 

XND 

XNDD  (1)  -XNDD  (30) 
X (1)  -X (30) 


Increment  (lbm/f t2/sec) 

Average  density  of  air  for  DX 
increment  (lbm/ft3) 

Viscosity  of  air  (lbf-sec/ f t2) 
Shear  force  over  the  increment  DX 
(lbf) 

Shear  force  over  the  condenser 
(lbf) 

Shear  force  over  the  evaporator 
and  adaiabatic 
region  (lbf) 

Summation  of  all  the  incremental 
shear  force 

values  over  the  entire  pipe  (lbf) 
Local  static  temperature  (R) 

Force  balance  check  (lbf) 

Wall  shear  stress  (lbf/ft2) 
Average  wall  shear  stress  for  the 
condenser 
(lbf/ft2) 

Total  air  temperature  (R) 

Average  axial  velocity  in 
condenser,  U  entering 
condenser  divided  by  2  (ft/sec) 
Average  velocity  at  beginning  of 
DX  increment 
(ft/sec) 

Average  velocity  at  end  of  DX 
increment  (ft/sec) 

Mass  flow  rate  at  the  upstream 
end  of  the  DX 
increment  (lbm/sec; 

Mass  flow  rate  at  the  downstream 
end  of  the  DX 
increment  (lbm/sec) 

Average  mass  flow  rate  for  DX 
increment  (lbm/sec) 

Mass  flow  rate  entering  condenser 
(lbm/sec) 

Local  axial  position  (ft) 
Non-dimensional  axial  position, 
XL/PL 

Non-dimensional  axial  position, 

X  (I)  /PL 

Axial  location  of  pressure 
measurements  (ft) 


**************************************************** 


*  INPUT  CONSTANTS  AND  DECLARE  TYPES  AND  DIMENSIONS  * 
**************************************************** 


IMPLICIT  REAL *8  (A-H,M,0-Z) 
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138 

REAL *8  X  (30)  ,P(30)  ,XNDD(30)  ,PND(30) 

139 

PL  = 

31.8125/12.0 

140 

X<1) 

=  0.0 

141 

X  (2) 

=  2.90625/12.0 

142 

X  (3) 

3  4.90625/12.0 

143 

X  (4) 

3  6.90625/12.0 

144 

X  (5) 

3  8.90625/12.0 

145 

X  (6) 

=  10.90625/12.0 

146 

X  (7) 

=  12.90625/12.0 

147 

X  (8) 

3  13.40625/12.0 

148 

X  (9) 

=  13.90625/12.0 

149 

XC10) 

=  14.40625/12.0 

150 

X(ll) 

3  14.90625/12.0 

151 

X  (12) 

=  15.53125/12.0 

152 

X  (1 3) 

=  16.28125/12.0 

153 

X  (14) 

=  16.90625/12.3 

154 

X  (15) 

=  17.40625/12.0 

155 

X  (16) 

=  17.90625/12.0 

156 

X  (17) 

=  18.90625/12.0 

157 

X  (1 8) 

=  19.90625/12.0 

158 

X  (1 9) 

=  20.90625/12.0 

159 

X  (20) 

=  21.90625/12.0 

160 

X  (21) 

=  22.90625/12.0 

161 

X  (22) 

=  23.90625/12.0 

162 

X  (2  3) 

3  24.90625/12.0 

163 

X  (24) 

=  25.90625/12.0 

164 

X  (25) 

=  26.90625/12.0 

165 

X  (26) 

3  27.90625/12.0 

166 

X  (27) 

=  28.90625/12.0 

167 

X  (28) 

=  29.90625/12.0 

168 

X  (29) 

=  30.90625/12.0 

169 

X  (30) 

=  31.8125/12.0 

170 

OPEN 

(UNIT=9 , STATUS= ' NEW' , FILE3 ' C : RESULTS . DAT 

171 

OPEN 

(UNIT=1 0 , STATUS= ' OLD’ , FILE= ' C : DATA . DAT ' ) 

172 

OPEN 

(UNIT=1 1 , STATUS= ’ NEW’ , FILE3 ' C :MACH . DAT '  ) 

173 

OPEN 

(UNIT=12 , STATUS3 ' NEW' .FILE3' C :REW .DAT’ ) 

174 

OPEN 

(UNIT= 1 3 , STATUS3 ’ NEW FILE 3 ’ C ; REX . DAT ’ ) 

175 

OPEN 

(UNIT=14 , STATUS3 ’ NEW’ .FILE3 ’ C : P2 . DAT ’ ) 

176 

OPEN 

(UNIT=15 , STATUS3 ' NEW' , FILE3' C :TW . DAT' ) 

177 

OPEN 

(UNIT3 1 6 , STATUS3 ’ NEW ' , FILE3 ' C : F . DAT ' ) 

178 

OPEN 

(UNIT3 17 , STATUS3 ' NEW’ , FILE3 ’ C : FBL . DAT ' ) 

179 

OPEN 

(UNIT=18, STATUS3' NEW' „ FILE3 ' C :RHOV2 . DAT ’ 

180 

OPEN 

(UNIT=1 9 , STATUS3 ’ NEW’ , FILE3 ' C : SF . DAT ' ) 

181 

OPEN 

(UNIT=20 , STATUS3 ’ NEW ’ , FILE3 ' C : DW . DAT ' ) 

182 

OPEN 

(UNIT=2 1 , STATUS3 ' NEW' , FILE3 ' C : PND . DAT ’ ) 

183 

OPEN 

(UNIT=22 .STATUS3 'NEW'  .FILE3' C :FREX .DAT’ ) 

184 

OPEN 

(UNIT=2  3 , STATUS3 ' NEW' , FILE3 ' C : FBLREX . DAT 

185 

OPEN 

(UNIT=24 .STATUS3' NEW’ .FILE3’ C :FBT .DAT' ) 

186 

OPEN 

(UNIT=25 , STATUS3 ' NEW FILE3 ’ C : FBTREX . DAT 

187 

OPEN 

(UNIT=26 , STATUS3 ' NEW' , FILE3 ' C : VERIFY . DAT 

188 

OPEN 

(UNIT=27 , STATUS3 ' NEW' , FILE3 ' C : FEVAP . DAT ' 

189 

READ 

(10,4000)  N 

B-5 


190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203  10 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216  C 

217  C 

218  C 

219  C 

220  C 

221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 


WRITE  (9,5020)  N 
READ  (10,4010)  PATM 
WRITE  (9,5030)  PATM 
READ  (10,4020)  TO 
WRITE  (9,5040)  TO 
READ  (10,4010)  PSINK 
WRITE  (9,5050)  PSINK 
DO  10  I  =  1,30 
READ  (10,4010)  P(I) 

P(I)  =  P (I)  *  144 
PND  (I)  =  P(I)/P(1) 

XNDD(I)  =  X  (I)  /PL 

WRITE  (21,5000)  XNDD(I)  ,PND(I) 

CONTINUE 
WRITE  (*,5060) 

READ  (*,4030)  FUD 
WRITE  (*,5061) 

READ  (*,4030)  MFFB 
WRITE  (*,5062) 

READ  (*,4030)  MFFS 

WRITE  (9,5070)  FUD,  MFFB,  MFFS 

PATM  =  PATM* 14 4 

PSINK  »  PSINK* 144 

DX  =  0.005 

A  =  3 . 592472958D8 

WRITE  (9,5080) 

********************************************* 

*  INITIALIZE  VARIABLES  AT  PIPE  END,  X  =  0.0  * 
********************************************* 


XL  =  0.0 
W  =0.0 
MIS  =  0.0 
Ml  =  0.0 
U1  =  0.0 
PI  =  P(l) 

P2  =  PI 
T  =  TO 
R  =  53.335 
GC  =  32.174 
PEX  =  PATM 
DWSUM  =0.0 
SFSUM  =0.0 
1  =  1 
PZ  =  P(I) 

PO  =  P(I  +  1) 
PT  =  P (1+2) 
XZ  =  X(I) 

XO  =  X(I+1) 
XT  =  X (1+2) 
XND  =0.0 
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242  CALL  SUBINC 

(FUD , A , TO , PEX , XND , P2 , DPS , RHOl , RHOV 1 , U1 , W1 , Ml , Ml S , 

24  3  + DWSUM , T , SFSUM , RH02 , RHOV2 , U2 , W2 , M2  S , M2 , DW , REX , REV , 

BET A, TV, SF, 

244  +F, FBL , FBT , PI , FREX , FBLREX , FBTREX , RHOBAR , UBAR , MFFB , 

MFFS) 

245  WRITE  (11,5000)  XND,  Ml 

246  WRITE  (12,5000)  XND,  REW 

247  WRITE  (13,5000)  XND,  Ml 

248  WRITE  (14,5000)  XND,  P2 

249  WRITE  (15,5000)  XND,  Ml 

250  WRITE  (16,5000)  XND,  Ml 

251  WRITE  (17,5000)  XND,  Ml 

252  WRITE  (18,5000)  XND,  RHOV1 

253  WRITE  (19,5000)  XND,  Ml 

254  WRITE  (20,5000)  XND,  DW 

255  WRITE  (22,5000)  XND,  Ml 

256  WRITE  (23,5000)  XND,  Ml 

257  C 

258  C  NOTE:  Ml  IS  USED  IN  ABOVE  WRITE  STATEMENTS  TO  ZERO 

VARIABLES  KNOWN 

259  C  TO  BE  ZERO  AT  END  OF  PIPE. 

260  C 

261  RHOV1  =  RHOV 2 

262  RHOl  =  RH02 

263  W1  =  W2 

264  PI  =  P2 

265  C 

266  C  ************************************ 

267  C  *****  BEGIN  MARCHING  DOWN  PIPE  ***** 

268  C  ************************************ 

269  C 

270  DO  200  J  =  1,81 

271  K  =  J 

272  C 

273  c  ************************************* 

274  c  ********  0.005  <  XL  <  0.405  ********* 

275  C  *****  0.001886  <  XND  <  0.152770  ***** 

276  C  ************************************* 

277  C 

278  XL  =  XL  +  DX 

279  CALL  SUBPRES  (PZ , PO , PT , XZ , XO, XT , XL , XND, P2) 

280  IF  (Ml  .LT.  0.01)  THEN 

281  CALL  SUBINC 

(FUD , A , TO , PEX , XND , P2 , DPS , RHOl , RHOV 1 , U1 , W1 , Ml , 

2  82  +M1 S , DWSUM , T , SFSUM , RH02 , RHOV 2 , U2 , W2 , M2  S , M2 , DW , REX , REW , 

283  +BETA , TW , SF , F , FBL , FBT , PI , FREX , FBLREX , FBTREX , RHOBAR , 

284  +UBAR, MFFB, MFFS) 

285  ELSE 

286  CALL  SUBCOMP 

(FUD, A, TO, PEX, XND, P2, DPS, PI ,RH01,RH0V1 ,U1,W1, 

287  +M1 , Ml S , DWSUM , T , SFSUM , RH02 , RHOV 2 , U2 , W2 , M2  S , M2 , DW , REX 
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,REW, 

288  +BETA ,TW , SF , F , FBL , FBT , FREX , FBLREX , FBTREX .RHOBAR , UBAR) 

289  ENDIF 

290  WRITE  (11,5000)  XND, M2 

291  WRITE  (12,5000)  XND.REW 

292  WRITE  (13,5000)  XND,REX 

293  WRITE  (14,5000)  XND.P2 

294  WRITE  (15,5000)  XND.TW 

295  WRITE  (16,5000)  XND, F 

296  WRITE  (17,5000)  XND, FBL 

297  WRITE  (18,5000)  XND.RH0V2 

298  WRITE  (19,5000)  XND.SF 

299  WRITE  (20,5000)  XND,DW 

300  WRITE  (22,5000)  XND, FREX 

301  WRITE  (23,5000)  XND, FBLREX 

302  IF  (XND  .GT.  0.099  .AND.  XND  . LT .  0.101)  THEN 

303  WRITE  (9,5011)  XND, F ,REW , REX , BET A, M2 , FREX , FBLREX 

304  ENDIF 

305  C 

306  C  ************************************************ 

307  C  *  CAPTURE  OF  VARIABLE  QUANTITIES  TO  BE  USED  TO  * 

308  C  *  VERIFIY  PROGRAM  BY  HAND  CALCULATION  * 

309  C  ************************************************ 

310  C 

311  IF  (K.EQ.3)  THEN 

312  PI  =  Pl/144.0 

313  P2  =  P2/144.0 


314 

RHOV1 

=  RHOV1/144.0 

315 

RHOV2 

=  RHOV2/144 . 0 

316 

RHOl  = 

RHOl/1728 . 0 

317 

RH02  = 

RH02/1728 , 0 

318 

RHOBAR 

=  RHOBAR/ 172 8.0 

319 

U1  =  Ul*12 . 0 

320 

U2  =  U2  *  1 2 . 0 

321 

UBAR  = 

UBAR* 12.0 

322 

DPS  = 

-DPS/20736.0 

323 

TW  =  TW/ 144.0 

324 

WRITE 

(26,*)  'K  =  3’ 

325 

WRITE 

(26,*)  ’Pi  =  '  ,P1 

326 

WRITE 

(26,*)  ’ P2  =  ' ,P2 

327 

WRITE 

(26,*)  ’ RHOV1  =  ' , 

RHOV1 

328 

WRITE 

(26,*)  'RHOV2  =  '  , 

RHOV2 

329 

WRITE 

(26,*)  'W1  =  ’ ,  W1 

330 

WRITE 

(26,*)  ’ W2  =  ' ,  W2 

331 

WRITE 

(26,*)  ’DW  =  ’ ,  DW 

332 

WRITE 

(26,*)  ' DWSUM  =  '  , 

DWSUM 

333 

WRITE 

(26,*)  'T  *  ’  ,  T 

334 

WRITE 

(26,*)  'RHOl  *  '  ,  ] 

RHOl 

335 

WRITE 

(26,*)  ’RH02  =  '  , 

RH02 

336 

WRITE 

(26,*)  ’RHOBAR  =  ' 

,  RHOBAR 

337 

WRITE 

(26,*)  'U1  =  ’ ,  U1 

338 

WRITE 

(26,*)  ' U2  =  ' ,  U2 

339 

WRITE  (26,*) 

'UBAR  = 

'  ,  UBAR 

340 

WRITE  (26,*) 

’Ml  =  ' 

,  Ml 

341 

WRITE  (26,*) 

'MIS  = 

'  ,  MIS 

342 

WRITE  (26,*) 

'M2  =  ’ 

,  M2 

34  3 

WRITE  (26,*) 

’  M2S  = 

'  ,  M2S 

344 

WRITE  (26,*) 

'ADPS  = 

’  ,  DPS 

345 

WRITE  (26,*) 

’  F  =  '  , 

F 

346 

WRITE  (26,*) 

'FREX  = 

' ,  FREX 

347 

WRITE  (26,*) 

'TW  =  ' 

,  TW 

348 

WRITE  (26,*) 

' SFSUM 

=  ' ,  SFSUM 

349 

WRITE  (26,*) 

'SF  =  ' 

,  SF 

350 

WRITE  (26,*) 

’  REW  = 

'  ,  REW 

351 

WRITE  (26,*) 

'REX  = 

'  ,  REX 

352 

PI  =  Pl*144 . 0 

353 

P2  =  P2*144 . 0 

354 

RHOV1  =  RH0V1*144 . 0 

355 

RHOV2  =  RHOV2* 144.0 

356 

RHOl  =  RH01*1728 . 0 

357 

RH02  =  RH02*1728 . 0 

358 

RHOBAR  =  RH0BAR*1728 

.0 

359 

U1  =  Ul/12.0 

360 

U2  =  U2/12.0 

361 

UBAR  =  UBAR/12 , 0 

362 

DPS  «  -DPS*20736 . 0 

363 

TW  =  TW* 144.0 

364 

END  IF 

365 

IF  (K.EQ.4)  THEN 

366 

PI  ■  Pl/144.0 

367 

P2  =  P2/144.0 

368 

RHOVl  =  RHOV1/144 . 0 

369 

RHOV2  =  RHOV2/ 144.0 

370 

RHOl  =  RHOl/1728 . 0 

371 

RH02  *  RH02/ 1728.0 

372 

RHOBAR  =  RHOBAR/ 172 8 

.0 

373 

U1  =  Ul *  12 . 0 

374 

U2  =  U2*12 . 0 

375 

UBAR  =  UBAR* 12.0 

376 

DPS  =  -DPS/20736.0 

377 

TW  =  TW/ 144.0 

378 

WRITE  (26,*) 

379 

WRITE  (26,*) 

’  K  =  4' 

380 

WRITE  (26,*) 

'PI  =  ’ 

,P1 

381 

WRITE  (26,*) 

’  P2  =  ' 

,P2 

382 

WRITE  (26,*) 

'RHOVl 

=  ' ,  RHOVl 

383 

WRITE  (26,*) 

'RHOV2 

=  '  ,  RHOV2 

384 

WRITE  (26,*) 

’  W1  =  ' 

,  W1 

385 

WRITE  (26,*) 

’  W2  =  ' 

,  W2 

386 

WRITE  (26,*) 

'DW  =  ' 

,  DW 

387 

WRITE  (26,*) 

’ DWSUM 

=  ' ,  DWSUM 

388 

WRITE  (26,*) 

’  T  =  ’  , 

T 

389 

WRITE  (26,*) 

'RHOl  * 

'  ,  RHOl 

390 

WRITE  (26,*) 

' RH02  = 

' ,  RH02 
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391 

WRITE  (26,*) 

' RHOBAR 

=  ' ,  RHOBAR 

392 

WRITE  (26,*) 

*U1  =  ' , 

U1 

393 

WRITE  (26,*) 

’  U2  =  '  , 

U2 

394 

WRITE  (26,*) 

'UBAR  = 

' ,  UBAR 

395 

WRITE  (26,*) 

'Ml  = 

Ml 

396 

WRITE  (26,*) 

'MIS  =  ’ 

,  MIS 

397 

WRITE  (26,*) 

'M2  =  '  , 

M2 

398 

WRITE  (26,*) 

’ M2S  =  ’ 

,  M2S 

399 

WRITE  (26,*) 

'ADPS  = 

'  .DPS 

400 

WRITE  (26,*) 

'F«'. 

F 

401 

WRITE  (26,*) 

'FREX  = 

* ,  FREX 

402 

WRITE  (26,*) 

’ TW  =  ’  , 

TW 

403 

WRITE  (26,*) 

' SFSUM  = 

' ,  SFSUM 

404 

WRITE  (26,*) 

’ SF  =  ' , 

SF 

405 

WRITE  (26,*) 

'REW  =  ’ 

,  REW 

406 

WRITE  (26,*) 

'REX  =  ' 

,  REX 

407 

PI  =  PI  *  1 44 . 0 

408 

P2  =  P2* 144.0 

409 

RHOV1  =  RHOV 1*144.0 

410 

RHOV2  =  RHOV2* 144.0 

411 

RHOl  =  RHOl* 1728.0 

412 

RH02  =  RH02* 1728.0 

413 

RHOBAR  =  RHOBAR* 1728 . 

0 

414 

U1  =  Ul/12.0 

415 

U2  =  U2/12.0 

416 

UBAR  =  UBAR/12 . 0 

417 

DPS  =  -DPS*20736 . 0 

418 

TW  =  TW* 144 . 0 

419 

END  IF 

420  C 

421  C 

*********************************** 

422  C 

*****  END  OF  VERIFICATION  ROUTINE  *' 

423  C 

*********************************** 

424  C 

425 

RHOV1  =  RHOV 2 

426 

RHOl  =  RH02 

427 

W1  =  W2 

428 

MIS  =  M2S 

429 

Ml  =  M2 

430 

U1  =  U2 

431 

PI  =  P2 

432  200 

CONTINUE 

433 

1  =  3 

4  34 

PZ  =  P(I) 

4  35 

PO  =  PCI+1) 

436 

PT  =  P(I+2) 

437 

XZ  =  X(I) 

438 

XO  =  X(I+1) 

439 

XT  =  X  (1+2) 

440 

DO  400  J  =  1,67 

441  C 
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442  C  ************************************* 

443  C  ********  0.410  <  XL  <  0.740  ********* 

444  C  *****  0.154656  <  XND  <  0.279136  ***** 

445  C  ************************************* 

446  C 

447  XL  =  XL  +  DX 

448  CALL  SUBPRES  (PZ , PO , PT , XZ , XO, XT , XL . XND, P2) 

449  IF  (Ml  .LT.  0.01)  THEN 

450  CALL  SUBINC 

(FUD , A ,T0 ,  PEX , XND , P2 , DPS ,RH01 .RHOV1 ,U1 ,W1 ,M1, 

451  +M1 S , DWSUM , T , SFSUM , RH02 , RHOV2 , U2 , W2 , M2  S , M2 , DW , REX , REW , 

452  +BETA , TW , SF , F , FBL , FBT , PI , FREX , FBLREX , FBTREX , RHOBAR , 

453  +UBAR , MFFB , MFFS) 

454  ELSE 

455  CALL  SUBCOMP 

(FUD , A , TO , PEX , XND , P2 , DPS , PI , RHOl ,RH0V1 ,U1.W1 , 

456  +M1 , MIS, DWSUM ,T. SFSUM ,RH02 .RHOV2 ,U2 , W2 ,M2S ,M2 , DW.REX 

.REW, 

457  +BETA , TW , SF , F , FBL , FBT , FREX , FBLREX , FBTREX , RHOBAR , UBAR) 

458  ENDIF 

459  WRITE  (11,5000)  XND, M2 

460  WRITE  (12,5000)  XND, REW 

461  WRITE  (13,5000)  XND, REX 

462  WRITE  (14,5000)  XND.P2 

463  WRITE  (15,5000)  XND.TW 

464  WRITE  (16,5000)  XND.F 

465  WRITE  (17,5000)  XND, FBL 

466  WRITE  (18,5000)  XND,RHOV2 

467  WRITE  (19,5000)  XND.SF 

468  WRITE  (20,5000)  XND.DW 

469  WRITE  (22,5000)  XND, FREX 

470  WRITE  (23,5000)  XND, FBLREX 

471  IF  (XND  ,GT .  0.199  .AND.  XND  .LT.  0.201)  THEN 

472  WRITE  (9,5011)  XND , F ,REW , REX , BETA ,M2 , FREX , FBLREX 

473  ENDIF 

474  RHOV1  =  RHOV2 

475  RHOl  =  RH02 

476  W1  =  W2 

477  MIS  »  M2S 

478  Ml  =  M2 

479  U1  =  U2 

480  PI  =  P2 

481  400  CONTINUE 

482  1  =  5 

483  PZ  =  P(I) 

484  PO  =  P(I+1) 

485  PT  =  P (1+2) 

486  XZ  =  X(I) 

487  XO  =  X(I+1) 

488  XT  =  X  (1+2) 

489  DO  600  J  =  1,67 

490  C 
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491  C  ************************************* 

492  C  ********  0.745  <  XL  <  1.075  ********* 

493  C  *****  0.281022  <  XND  <  0.405501  ***** 

494  C  ************************************* 

495  C 

496  XL  =  XL  +  DX 

497  CALL  SUBPRES  (PZ , PO , PT , XZ , XO , XT , XL , XND, P2) 

498  IF  (Ml  .LT.  0.01)  THEN 

499  CALL  SUBINC 

(FUD , A ,T0 , PEX , XND , P2 , DPS ,RH01 .RHOV1 ,U1 ,V1 ,M1 , 

500  +M1 S , DWSUM , T , SFSUM , RH02 , RHOV2 , U2 , W2 , M2  S , M2 , DW , REX , REV , 

501  +BETA , TV , SF , F , FBL , FBT , PI , FREX , FBLREX , FBTREX , RHOBAR , 

502  +UBAR , MFFB , MFFS) 

503  ELSE 

504  CALL  SUBCOMP 

(FUD , A .TO , PEX , XND , P2 , DPS , PI ,RH01 ,RH0V1 ,U1 ,W1 . 

505  +M1 ,M1S,DVSUM,T,SFSUM,RH02,RHOV2 ,U2,V2,M2S,M2 ,DV,REX 

.REV, 

506  +BETA , TV , SF , F , FBL , FBT , FREX , FBLREX , FBTREX , RHOBAR , UBAR) 

507  ENDIF 


508 

VRITE  (11,5000) 

XND, M2 

509 

!  VRITE  (12,5000) 

XND, REV 

510 

VRITE  (13,5000) 

XND, REX 

511 

VRITE  (14,5000) 

XND, P2 

512 

VRITE  (15,5000) 

XND, TV 

513 

VRITE  (16,5000) 

XND ,  F 

514 

VRITE  (17,5000) 

XND, FBL 

515 

VRITE  (18,5000) 

XND,RHOV2 

516 

VRITE  (19,5000) 

XND, SF 

517 

VRITE  (20,5000) 

XND.DV 

518 

VRITE  (22,5000) 

XND, FREX 

519 

VRITE  (23,5000) 

XND, FBLREX 

520 

IF  (XND  .GT.  0.298  .AND.  XND  .LT.  0.301)  THEN 

521 

VRITE  (9,5011) 

XND , F , REV , REX , BET A , M2 , FREX , FBLREX 

522 

ENDIF 

523 

IF  (XND  .GT.  0. 

398  .AND.  XND  .LT.  0.401)  THEN 

524 

VRITE  (9,5011) 

XND , F , REV , REX , BETA , M2 , FREX , FBLREX 

525 

ENDIF 

526 

RHOV1  =  RHOV2 

527 

RHOl  =  RH02 

528 

VI  =  V  2 

529 

MIS  =  M2S 

530 

Ml  =  M2 

531 

U1  =  U2 

532 

PI  =  P2 

533  600 

CONTINUE 

534 

1  =  7 

535 

PZ  =  P(I) 

536 

PO  =  P(I+1) 

537 

PT  =  P  (1+2) 

538 

XZ  =  X(I) 

539 

XO  =  X(I+1) 
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540 

541 

542  C 

543  C 

544  C 

545  C 

546  C 

547  C 

548 

549 

550 

551 

552 

553 

554 

555 

556 

557 

558 

559 

560 

561 

562 

563 

564 

565 

566 

567 

568 

569 

570 

571 

572 

573 

574 

575 

576 

577 

578 

579  800 

580 

581 

582 

583 

584 

585 

586 

587 

588  C 


XT  =  X  (1+2) 

DO  800  J  =  1,16 

************************************* 

********  1.080  <  XL  <  1.155  ********* 

*****  0.407387  <  XND  <  0.435678  ***** 
************************************* 

XL  =  XL  +  DX 

CALL  SUBPRES  (PZ , PO , PT , XZ , XO , XT , XL , XND , P2) 

IF  (Ml  .LT.  0.01)  THEN 
CALL  SUBINC 

(FUD , A , TO , PEX , XND , P2 , DPS , RH01 ,RH0V1 ,U1 ,W1 ,M1 , 

+M1S , DWSUM.T , SFSUM , RH02 , RHOV  2 , U2 , W2 , M2  S , M2 , DW , REX , REW , 
+BETA , TW , SF , F , FBL , FBT , Pi , FREX , FBLREX , FBTREX , RHOBAR , 
+UBAR ,MFFB,MFFS) 

ELSE 

CALL  SUBCOMP 

(FUD , A , TO , PEX , XND , P2 , DPS , PI , RH01 ,RHOV 1 , U1 , W1 , 

+M1, MIS, DWSUM.T, SFSUM, RH02 , RH0V2 , U2 , W2 ,M2S , M2, DW, REX 
,REW , 

+BETA , TW , SF , F , FBL , FBT , FREX , FBLREX , FBTREX , RHOBAR , UBAR) 
ENDIF 

WRITE  (11,5000)  XND, M2 
WRITE  (12,5000)  XND, REW 
WRITE  (13,5000)  XND, REX 
WRITE  (14,5000)  XND,P2 
WRITE  (15,5000)  XND.TW 
WRITE  (16,5000)  XND.F 
WRITE  (17,5000)  XND, FBL 
WRITE  (18,5000)  XND,RHOV2 
WRITE  (19,5000)  XND,SF 
WRITE  (20,5000)  XND.DW 
WRITE  (22,5000)  XND, FREX 
WRITE  (23,5000)  XND, FBLREX 
RHOV1  =  RHOV 2 
RHOl  =  RH02 
W1  =  W2 
MIS  =  M2S 
Ml  =  M2 
U1  =  U2 
PI  =  P2 
CONTINUE 
1  =  9 
PZ  =  P(I) 

PO  *  P(I+1) 

PT  =  P(I+2) 

XZ  =  X(I) 

XO  =  X(I  +  1) 

XT  =  X  (1+2) 

DO  1000  J  =  1,17 
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589  C  ************************************* 

590  C  ********  1.16O  <  XL  <  1.240  ********* 

591  C  *****  0.437564  <  XND  <  0.467741  ***** 

592  C  ************************************* 

^93  C 

594  XL  *  XL  +  DX 

595  CALL  SUBPRES  (PZ , PO, PT , XZ , XO, XT , XL , XND, P2) 

596  IF  (Ml  .LT.  0.01)  THEN 

597  CALL  SUBINC 

(FUD , A , TO , PEX , XND ,  P2 , DPS , RHOl , RHOV 1 ,U1 ,W1 ,M1 , 

598  +M1S , DWSUM,T , SFSUM ,RH02 , RHOV2 , U2 , W2 ,M2S ,M2 , DW ,REX ,REW , 

599  +BETA , TW , SF , F , FBL , FBT , PI , FREX , FBLREX , FBTREX , RHOBAR , 

600  +UBAR , MFFB , MFFS) 

601  ELSE 

602  CALL  SUBCOMP 

(FUD, A, TO, PEX, XND , P2 , DPS, PI ,RH01 ,RH0V1 ,U1 ,W1 , 

603  +M1 , Ml S , DWSUM , T , SFSUM , RH02 , RH0V2 , U2 , W2 , M2  S , M2 , DW , REX 

,REW, 

604  +BETA , TW , SF , F , FBL , FBT , FREX , FBLREX , FBTREX , RHOBAR , UBAR) 

605  END IF 

606  WRITE  (11,5000)  XND, M2 

607  WRITE  (12,5000)  XND,REW 

608  WRITE  (13,5000)  XND, REX 

609  WRITE  (14,5000)  XND.P2 

610  WRITE  (15,5000)  XND.TW 

611  WRITE  (16,5000)  XND, F 

612  WRITE  (17,5000)  XND, FBL 

613  WRITE  (18,5000)  XND, RHOV 2 

614  WRITE  (19,5000)  XND,SF 

615  WRITE  (20,5000)  XND.DW 

616  WRITE  (22,5000)  XND, FREX 

617  WRITE  (23,5000)  XND, FBLREX 

618  RHOV1  =  RHOV2 

619  RHOl  =  RH02 

620  W1  =  W2 

621  MIS  =  M2S 

622  Ml  *  M2 

623  U1  *  U2 

624  PI  =  P2 

625  1000  CONTINUE 

626  I  =  11 

627  PZ  =  P(I) 

628  PO  *  P(I+1) 

629  PT  =  P(I+2) 

630  XZ  »  X(I) 

631  XO  =  X(I+1) 

632  XT  =  X  (1+2) 

633  DO  1200  J  =  1,23 

634  C 

635  C  ************************************* 

636  c  ********  1.245  <  XL  <  1.355  ********* 

637  C  *****  0.469627  <  XND  <  0.511120  ***** 
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I************************************ 


638  C 

639  C 

640  XL  =  XL  +  DX 

641  CALL  SUBPRES  (PZ , PO , PT , XZ , XO, XT , XL , XND. P2) 

642  IF  (XND  .GT.  0.491  .AND.  XND  .LT.  0.508)  PEX  =  P2 

643  IF  (XND  .GT.  0.508)  PEX  =  PSINK 

644  IF  (Ml  .LT.  0.01)  THEN 

645  CALL  SUBINC 

(FUD , A , TO , PEX . XND , P2 , DPS , RH01 , RHOV1 , U1 , W1 , Ml , 

64  6  +M1S ,  DWSlfM ,  T ,  SFSUM ,  RH02  ,  RHOV2  ,  U2  ,  W2  ,  M2  S ,  M2  ,  DW ,  REX  ,  REW , 

647  +  BETA , TW , SF , F , FBL , FBT , P 1 , FREX , FBLREX , FBTREX , RHOBAR , 

648  +UBAR , MFFB , MFFS) 

649  ELSE 

650  CALL  SUBCOMP 

(FUD , A , TO , PEX , XND , P2 , DPS , PI , RHOl ,RHOV 1 , U1 , W1 , 

651  +M1 , Ml S , DWSUM , T , SFSUM , RH02 , RHOV2 , U2 , W2 , M2  S , M2 , DW , REX 

,REV, 

652  +BETA , TW , SF , F , FBL , FBT , FREX , FBLREX , FBTREX , RHOBAR . UBAR) 

653  ENDIF 

654  WRITE  (11,5000)  XND, M2 

655  WRITE  (12,5000)  XND, REW 

656  WRITE  (13,5000)  XND, REX 

657  WRITE  (14,5000)  XND.P2 

658  WRITE  (15,5000)  XND, TV 

659  WRITE  (16,5000)  XND, F 

660  WRITE  (17,5000)  XND, FBL 

661  WRITE  (18,5000)  XND,RHOV2 

662  WRITE  (19,5000)  XND,SF 

663  WRITE  (20,5000)  XND,DW 

664  WRITE  (22,5000)  XND, FREX 

665  WRITE  (23,5000)  XND, FBLREX 

666  IF  (XND  .GT.  0.508)  THEN 

667  WRITE  (24,5000)  XND, FBT 

668  WRITE  (25,5000)  XND, FBTREX 

669  ENDIF 

670  IF  (XND  .GT.  0.506  .AND.  XND  .LT.  0.508)  THEN 

671  UBARC  =  UBAR/ 2 

672  REXBRC  =  REX/ 2 

673  ENDIF 

674  IF  (XND  .GT.  0.508  .AND.  XND  .LT.  0.510)  THEN 

675  SFEA  =  SFSUM 

676  RHOC  =  RH02 

677  ENDIF 

678  RHOV1  =  RHOV2 

679  RHOl  *  RH02 

680  W1  =  W2 

681  MIS  =  M2S 

682  Ml  =  M2 

683  U1  =  U2 

684  PI  *  P2 

685  IF  (XND  .GT.  0.489  .AND.  XND  .LT.  0.491)  WCK  =  DWSUM 

686  1200  CONTINUE 
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687  I  =  13 

688  PZ  =  PCI) 

689  PO  =  PCI+l) 

690  PT  =  P (1+2) 

691  XZ  =  XCI) 

692  XO  =  XCI  +  l) 

693  XT  =  X  (1+2) 

694  DO  1400  J  =  1,19 

695  C 

696  C  ************************************* 

697  C  ********  1.360  <  XL  <  1.450  ********* 

698  C  *****  0.513006  <  XND  <  0.546955  ***** 

699  C  ************************************* 

700  C 

701  XL  =  XL  +  DX 

702  CALL  SUBPRES  (PZ , PO , PT , XZ , XO , XT , XL , XND , P2) 

703  IF  (Ml  .LT.  0.01)  THEN 

704  CALL  SUBINC 

(FUD , A , TO , PEX , XND , P2 , DPS , RH01 ,RH0V1 ,U1 ,W1 ,M1, 

705  +M1S , DWSUM , T , SFSUM , RH02 , RH0V2 , U2 , W2 , M2  S , M2 , DW , REX , REW , 

706  +BETA , TW , SF , F , FBL , FBT , PI , FREX , FBLREX , FBTREX , RHOBAR , 

707  +UBAR ,MFFB,MFFS) 

708  ELSE 

709  CALL  SUBCOMP 

(FUD , A, TO, PEX, XND, P2, DPS, PI ,RH01 ,RH0V1 ,U1 ,W1 , 

710  +M1 , MIS, DWSUM, T, SFSUM, RH02 ,RHOV2 , U2 ,W2 ,M2S ,M2 , DW ,REX 

,REW, 

711  +BETA , TW , SF , F , FBL , FBT , FREX , FBLREX , FBTREX , RHOBAR , UBAR) 

712  ENDIF 

713  WRITE  (11,5000)  XND, M2 

714  WRITE  (12,5000)  XND, REW 

715  WRITE  (13,5000)  XND, REX 

716  WRITE  (14,5000)  XND.P2 

717  WRITE  (15,5000)  XND.TW 

718  WRITE  (16,5000)  XND,F 

719  WRITE  (17,5000)  XND, FBL 

720  WRITE  (18,5000)  XND.RHOV2 

721  WRITE  (19,5000)  XND,SF 

722  WRITE  (20,5000)  XND.DW 

723  WRITE  (22,5000)  XND, FREX 

724  WRITE  (23,5000)  XND, FBLREX 

725  WRITE  (24,5000)  XND, FBT 

726  WRITE  (25,5000)  XND, FBTREX 

727  RHOV1  =  RHOV2 

728  RHOl  =  RH02 

729  W1  »  W2 

730  MIS  =  M2S 

731  Ml  =  M2 

732  U1  =  U2 

733  PI  =  P2 

734  1400  CONTINUE 

735  I  =  15 
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736  PZ  =  PCI) 

737  PO  =  PCI  +  1) 

738  PT  =  P (1+2) 

739  XZ  =  X(I) 

740  XO  =  X(I+1> 

741  XT  =  X  (1+2) 

742  DO  1600  J  =  1,25 

743  C 

744  c  ************************************* 

745  C  ********  1.455  <  XL  <  1.575  ********* 

746  C  *****  0.548841  <  XND  <  0.594106  ***** 

747  C  ************************************* 

748  C 

749  XL  =  XL  +  DX 

750  CALL  SUBPRES  (PZ , PO , PT , XZ , XO , XT , XL , XND , P2) 

751  IF  (Ml  .LT.  0.01)  THEN 

752  CALL  SUBINC 

(FUD , A , TO , PEX , XND , P2 , DPS ,RH01 .RH0V1 ,U1 ,W1 ,M1 , 

753  +M1S , DWSUM , T , SFSUM ,RH02 ,RH0V2 , U2 , W2 ,M2S ,M2 , DW ,REX ,REW , 

754  +BETA , TW , SF , F , FBL , FBT , PI , FREX , FBLREX , FBTREX , RHOBAR . 

755  +UBAR , MFFB , MFFS) 

756  ELSE 

757  CALL  SUBCOMP 


(FUD , A , TO , PEX , XND , P2 , DPS , PI , RHOl , RHOV 1 , U1 , W1 , 

758  +M1 , Ml S , DWSUM , T , SFSUM , RH02 , RHOV2 , U2 , W2 , M2  S , M2 , DW , REX 

,REW , 

75  9  +  BETA , TW , SF , F , FBL , FBT , FREX , FBLREX , FBTREX , RHOBAR , UBAR) 


760 

ENDIF 

761 

WRITE 

(11,5000) 

XND, M2 

762 

WRITE 

(12,5000) 

XND.REW 

763 

WRITE 

(13,5000) 

XND, REX 

764 

WRITE 

(14,5000) 

XND , P2 

765 

WRITE 

(15,5000) 

XND , TW 

766 

WRITE 

(16,5000) 

XND,  F 

767 

WRITE 

(17,5000) 

XND, FBL 

768 

WRITE 

(18,5000) 

XND , RHOV2 

769 

WRITE 

(19,5000) 

XND , SF 

770 

WRITE 

(20,5000) 

XND , DW 

771 

WRITE 

(22,5000) 

XND , FREX 

772 

WRITE 

(23,5000) 

XND, FBLREX 

773 

WRITE 

(24,5000) 

XND, FBT 

774 

WRITE 

(25,5000) 

XND, FBTREX 

775 

RHOV1 

=  RHOV 2 

776 

RHOl  = 

RH02 

777 

W1  =  W2 

778 

MIS  = 

M2S 

779 

Ml  =  M2 

780 

U1  =  U2 

781 

PI  =  P2 

782  1600 

CONTINUE 

783 

I  =  17 

784 

PZ  =  P(I) 

B-17 


785 

786 

787 

788 

789 

790 

791  C 

792  C 

793  C 

794  C 

795  C 

796  C 

797 

798 

799 

800 

801 

802 

803 

804 

805 

806 

807 

808 

809 

810 
811 
812 

813 

814 

815 

816 

817 

818 

819 

820 
821 
822 
82  3 

824 

825 

826 

827 

828 

829 

830 

831 

832 


PO  =  PCI+1) 

PT  =  P  (1+2) 

XZ  =  X(I) 

XO  =  XCI  +  1) 

XT  =  X<I+2) 

DO  1800  J  =  1,33 

********  1.580  <  XL  <  1.740  ********* 
*****  0.595992  <  XND  <  0.656346  ***** 


XL  =  XL  +  DX 

CALL  SUBPRES  (PZ , PO , PT , XZ , XO , XT , XL , XND , P2) 

IF  (Ml  .LT.  0.01)  THEN 
CALL  SUBINC 

(FUD , A ,T0 , PEX , XND , P2 , DPS ,RH01 ,RH0V1 ,U1 ,W1 ,M1 , 

+M1 S , DWSUM , T , SFSUM , RH02 , RHOV2 , U2 , W2 , M2  S , M2 , DW , REX , REW , 
+BETA , TW , SF , F , FBL , FBT , PI , FREX , FBLREX , FBTREX , RHOBAR , 
+UBAR,MFFB,MFFS) 

ELSE 

CALL  SUBCOMP 


(FUD . A , TO , PEX , XND , P2 , DPS , PI , RHOl , RHOV 1 , U1 , W1 , 

+M1 , Ml S , DWSUM , T , SFSUM , RH02 , RHOV 2 , U2 , W2 , M2  S , M2 , DW , REX 
,REW, 

♦BETA , TW , SF , F , FBL , FBT , FREX , FBLREX , FBTREX , RHOBAR , UBAR) 
ENDIF 

WRITE  (11,5000)  XND, M2 
WRITE  (12,5000)  XND, REW 
WRITE  (13,5000)  XND, REX 
WRITE  (14,5000)  XND,P2 
WRITE  (15,5000)  XND,TW 
WRITE  (16,5000)  XND,F 
WRITE  (17,5000)  XND, FBL 
WRITE  (18,5000)  XND.RHOV2 
WRITE  (19,5000)  XND,SF 
WRITE  (20,5000)  XND.DW 
WRITE  (22,5000)  XND, FREX 
WRITE  (23,5000)  XND, FBLREX 
WRITE  (24,5000)  XND, FBT 
WRITE  (25,5000)  XND, FBTREX 

IF  (XND  .GT.  0.599  .AND.  XND  .LT.  0.601)  THEN 
WRITE  (9,5010) 

XND , F , REW , REX , BETA , M2 , FREX , FBLREX , FBTREX 
ENDIF 

RHOV1  =  RHOV2 
RHOl  =  RH02 
W1  =  W2 
MIS  =  M2S 
Ml  =  M2 
U1  =  U2 
PI  =  P2 


B-18 


833  1800  CONTINUE 

834  I  =  19 


835 

PZ 

P(I) 

836 

PO 

= 

P(I  +  1) 

837 

PT 

St 

P  (1+2) 

838 

XZ 

= 

X(I) 

839 

XO 

X(I  +  1) 

840 

XT 

= 

X  (1+2) 

841 

DO 

2000  J  = 

842  C 

843  C  ************************************* 

844  C  ********  1.745  <  XL  <  1.905  ********* 

845  C  *****  0.658232  <  XND  <  0.718585  ***** 

846  C  ************************************* 

847  C 

848  XL  =  XL  +  DX 

849  CALL  SUBPRES  (PZ , PO , PT , XZ , XO . XT , XL , XND , P2) 

850  IF  (Ml  .LT.  0.01)  THEN 

851  CALL  SUBINC 

(FUD , A , TO , PEX , XND , P2 , DPS , RHOl , RHOV 1 , U1 , W1 , Ml , 

852  +M1S , DWSUM , T , SFSUM , RH02 , RHOV 2 , U2 , W2 , M2  S , M2 . DW , REX , REW , 

85  3  +BETA , TW , SF , F , FBL , FBT „ PI , FREX , FBLREX , FBTREX , RHOBAR , 

854  +UBAR,MFFB,MFFS) 

855  ELSE 

856  CALL  SUBCOMP 

(FUD , A , TO , PEX , XND , P2 , DPS , PI , RHOl ,RHOV 1 , U1 , W1 , 

857  +M1 , Ml S , DWSUM , T , SFSUM , RH02 , RHOV 2 , U2 , W2 , M2  S , M2 , DW , REX 

.REW, 

85  8  +BETA , TW , SF , F , FBL , FBT , FREX , FBLREX , FBTREX , RHOBAR . UBAR) 

859  ENDIF 

860  WRITE  (11,5000)  XND, M2 

861  WRITE  (12,5000)  XND, REW 

862  WRITE  (13,5000)  XND, REX 

863  WRITE  (14,5000)  XND.P2 

864  WRITE  (15,5000)  XND.TW 

865  WRITE  (16,5000)  XND, F 

866  WRITE  (17,5000)  XND, FBL 

867  WRITE  (18,5000)  XND.RHOV2 

868  WRITE  (19,5000)  XND.SF 

869  WRITE  (20,5000)  XND.DW 

870  WRITE  (22,5000)  XND, FREX 

871  WRITE  (23,5000)  XND, FBLREX 

872  WRITE  (24,5000)  XND, FBT 

873  WRITE  (25,5000)  XND, FBTREX 

874  IF  (XND  .GT.  0.698  .AND.  XND  .LT.  0.701)  THEN 

875  WRITE  (9,5010) 

XND , F , REW , REX , BETA , M2 , FREX , FBLREX , FBTREX 

876  ENDIF 

877  RHOV1  =  RHOV 2 

878  RHOl  =  RH02 

879  W1  =  W2 

880  MIS  *  M2S 
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881  Ml  =  M2 

882  U1  =  U2 

883  PI  =  P2 

884  2000  CONTINUE 

885  I  =  21 

886  PZ  =  PCI) 

887  PO  =  P(I  +  1) 

888  PT  =  P (1+2) 

889  XZ  =  X(I) 

890  XO  -  XCI+1) 

891  XT  =  X  (1+2) 

892  DO  2200  J  =  1,34 

893  C 

894  c  ************************************* 

895  C  ********  1.910  <  XL  <  2.075  ********* 

896  C  *****  0.720472  <  XND  <  0.782711  ***** 

897  C  ************************************* 

898  C 

899  XL  =  XL  +  DX 

900  CALL  SUBPRES  (PZ , PO, PT , XZ , XO, XT , XL , XND, P2) 

901  IF  (Ml  .LT.  0.01)  THEN 

902  CALL  SUBINC  , 

(FUD , A , TO , PEX , XND , P2 , DPS , RHOl , RHOV1 , U1 , W1 , Ml , 

903  +M1S , DWSUM , T , SFSUM , RH02 , RHOV2 , U2 , W2 , M2  S , M2 , DW , REX , REW , 

904  +BETA , TW , SF , F , FBL , FBT , PI , FREX , FBLREX , FBTREX , RHOBAR , 

905  +UBAR , MFFB , MFFS) 

906  ELSE 

907  CALL  SUBCOMP 

(FUD, A, TO, PEX , XND, P2 , DPS, PI ,RH01 ,RH0V1 ,U1 ,W1 , 

908  +M1, MIS, DWSUM, T, SFSUM, RH02 ,RHOV2 ,U2 , W2 ,M2S , M2, DW, REX 

,REW , 

909  +BETA , TW , SF , F , FBL , FBT , FREX , FBLREX . FBTREX .RHOBAR , UBAR) 

910  END IF 

911  WRITE  (11,5000)  XND, M2 

912  WRITE  (12,5000)  XND, REW 

913  WRITE  (13,5000)  XND, REX 

914  WRITE  (14,5000)  XND,P2 

915  WRITE  (15,5000)  XND,TW 

916  WRITE  (16,5000)  XND, F 

917  WRITE  (17,5000)  XND, FBL 

918  WRITE  (18,5000)  XND.RHOV2 

919  WRITE  (19,5000)  XND.SF 

920  WRITE  (20,5000)  XND.DW 

921  WRITE  (22,5000)  XND, FREX 

922  WRITE  (23,5000)  XND, FBLREX 

923  WRITE  (24,5000)  XND, FBT 

924  WRITE  (25,5000)  XND, FBTREX 

925  RHOV1  =  RHOV2 

926  RHOl  *  RH02 

927  W1  *  W2 

928  MIS  =  M2S 

929  Ml  =  M2 
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930  U1  =  U2 

931  PI  =  P2 

932  2200  CONTINUE 

933  I  =  23 

934  PZ  =  PCI) 

935  PO  =  P(I+1) 

936  PT  =  P (1+2) 

937  XZ  =  X(I) 

938  XO  =  X(I+1> 

939  XT  =  X (1+2) 

940  DO  2400  J  =  1,33 

941  C 

942  C  *********************************** 

943  C  ********  2.080  <  XL  <  2.240  ******* 

944  c  *****  0.784597  <  XND  <  844951  ***** 

945  C  *********************************** 

946  C 

947  XL  =  XL  +  DX 

948  CALL  SUBPRES  (PZ , PO , PT , XZ , XO, XT , XL , XND , P2) 

949  IF  (Ml  .LT.  0.01)  THEN 

950  CALL  SUBINC 

(FUD , A , TO , PEX , XND , P2 , DPS , RHOl , RHOV 1 , U1 , W1 , Ml , 

951  +M1S , DWSUM , T , SFSUM , RH02 , RHOV 2 , U2 , W2 , M2  S , M2 , DW , REX , REW , 

952  +BETA , TW , SF , F , FBL , FBT , PI , FREX , FBLREX , FBTREX ,RHOBAR , 

953  +UBAR , MFFB , MFFS) 

954  ESE 

955  CALL  SUBCOMP 

(FUD , A , TO , PEX , XND , P2 , DPS , PI , RHOl , RHOV 1 , U1 , W1 , 

956  +M1 ,M1S , DWSUM , T , SFSUM , RH02 , RHOV2 , U2 , W2 , M2  S , M2 , DW , REX 

.REW, 

957  +BETA , TW , SF , F , FBL , FBT , FREX , FBLREX , FBTREX , RHOBAR , UBAR) 

958  ENDIF 

959  WRITE  (11,5000)  XND, M2 

960  WRITE  (12,5000)  XND, REW 

961  WRITE  (13,5000)  XND, REX 

962  WRITE  (14,5000)  XND.P2 

963  WRITE  (15,5000)  XND,TW 

964  WRITE  (16,5000)  XND, F 

965  WRITE  (17,5000)  XND, FBL 

966  WRITE  (18,5000)  XND.RHOV2 

967  WRITE  (19,5000)  XND.SF 

968  WRITE  (20,5000)  XND.DW 

969  WRITE  (22,5000)  XND, FREX 

970  WRITE  (23,5000)  XND, FBLREX 

971  WRITE  (24,5000)  XND, FBT 

972  WRITE  (25,5000)  XND, FBTREX 

973  IF  (XND  .GT.  0.799  .AND.  XND  .LT.  0.801)  THEN 

974  WRITE  (9,5010) 

XND , F , REW , REX , BETA , M2 , FREX , FBLREX , FBTREX 

975  ENDIF 

976  RHOV1  =  RHOV 2 

977  RHOl  =  RH02 
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978  W1  =  W2 

979  MIS  =  M2S 

980  Ml  =  M2 

981  U1  =  U2 

982  PI  =  P2 

983  2400  CONTINUE 

984  I  =  25 

985  PZ  =  P(I) 

986  PO  =  PCI  +  l) 

987  PT  =  P (1+2) 

988  XZ  =  X(I) 

989  XO  =  XCI+1) 

990  XT  =  X  (1+2) 

991  DO  2600  J  =  1,33 

992  C 

993  C  ************************************* 

994  C  ********  2.245  <  XL  <  2.405  ********* 

995  C  *****  0.846837  <  XND  <  0.907191  ***** 

996  C  ************************************* 

997  C 

998  XL  =  XL  +  DX 

999  CALL  SUBPRES  (PZ , PO, PT , XZ , XO , XT , XL , XND, P2) 

1000  IF  (Ml  .LT.  0.01)  THEN 

1001  CALL  SUBINC 

(FUD , A , TO , PEX , XND, P2 , DPS ,RH01 ,RH0V1 , U1 , W1 ,M1 , 

1002  +M1S , DWSUM , T , SFSUM , RH02 , RH0V2 , U2 , W2 , M2  S , M2 , DW , REX , REW , 

1003  + BET A , TW , SF , F , FBL , FBT , P 1 , FREX , FBLREX , FBTREX , RHOBAR , 

1004  +UBAR,MFFB,MFFS) 

1005  ELSE 

1006  CALL  SUBCOMP 

(FUD, A, TO, PEX, XND, P2, DPS, PI ,RH01 .RHOV1 ,U1 ,W1 , 

1007  +M1 , Ml S , DWSUM , T , SFSUM , RH02 , RHOV2 , U2 , W2 , M2  S , M2 , DW , REX 

.REW, 

1 008  +BETA , TW , SF , F , FBL , FBT , FREX , FBLREX , FBTREX , RHOBAR , UBAR) 

1009  ENDIF 

1010  WRITE  (11,5000)  XND, M2 

1011  WRITE  (12,5000)  XND, REW 

1012  WRITE  (13,5000)  XND, REX 

1013  WRITE  (14,5000)  XND,P2 

1014  WRITE  (15,5000)  XND.TW 

1015  WRITE  (16,5000)  XND, F 

1016  WRITE  (17,5000)  XND, FBL 

1017  WRITE  (18,5000)  XND.RHOV2 

1018  WRITE  (19,5000)  XND.SF 

1019  WRITE  (20,5000)  XND.DW 

1020  WRITE  (22,5000)  XND, FREX 

1021  WRITE  (23,5000)  XND, FBLREX 

1022  WRITE  (24,5000)  XND, FBT 

1023  WRITE  (25,5000)  XND, FBTREX 

1024  IF  (XND  .GT.  0.899  .AND.  XND  .LT.  0.901)  THEN 

1025  WRITE  (9,5010) 

XND , F , REW , REX , BETA , M2 , FREX , FBLREX , FBTREX 
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1026 

1027 

1028 

1029 

1030 

1031 

1032 

1033 

1034 

1035 

1036 

1037 

1038 

1039 

1040 

1041 

1042 

1043 

1044 

1045 

1046 

1047 

1048 

1049 

1050 

1051 

1052 

1053 

1054 

1055 

1056 

1057 

1058 

1059 

1060 
1061 
1062 

1063 

1064 

1065 

1066 

1067 

1068 

1069 

1070 

1071 

1072 

1073 


ENDIF 

RHOV1  =  RHOV2 
RHOl  =  RH02 
W1  =  W2 
MIS  =  M2S 
Ml  =  M2 
U1  =  U2 
PI  =  P2 
2600  CONTINUE 
I  =  27 
PZ  =  PCI) 

PO  =  P(I  +  1) 

PT  =  P (1+2) 

XZ  =  X(I) 

XO  =  X(I  +  1) 

XT  =  X  (1+2) 

DO  2800  J  =  1,34 
C 

C  ********  2.410  <  XL  <  2.575  ********* 

C  *****  0.909077  <  XND  <  0.971316  ***** 

C 


XL  =  XL  +  DX 

CALL  SUBPRES  (PZ , PO , PT , XZ , XO , XT , XL , XND, P2) 

IF  (Ml  .LT.  0.01)  THEN 
CALL  SUBINC 

(FUD , A , TO , PEX . XND , P2 , DPS , RHOl , RHOV 1 , U1 , W1 , Ml , 

+M1 S , DWSUM , T , SFSUM , RH02 , RHOV 2 , U2 , W2 , M2  S , M2 , DW , REX 
,REW, 

+BETA , TW , SF , F , FBL , FBT , PI , FREX , FBLREX , FBTREX , RHOBAR , 
+UBAR , MFFB , MFFS) 

ELSE 

CALL  SUBCOMP 

(FUD, A, TO, PEX, XND, P2, DPS, PI, RHOl ,RH0V1 , U1 , W1 , 

+M1 , Ml S , DWSUM , T , SFSUM , RH02 , RHOV 2 , U2 , W2 , M2  S , M2 , DW , 
REX.REW, 

+BETA , TW , SF , F , FBL , FBT , FREX , FBLREX , FBTREX , RHOBAR , UBAR) 


ENDIF 

WRITE  (11,5000) 
WRITE  (12,5000) 
WRITE  (13,5000) 
WRITE  (14,5000) 
WRITE  (15,5000) 
WRITE  (16,5000) 
WRITE  (17,5000) 
WRITE  (18,5000) 
WRITE  (19,5000) 
WRITE  (20,5000) 
WRITE  (22,5000) 
WRITE  (23,5000) 
WRITE  (24,5000) 


XND, M2 
XND,REW 
XND, REX 
XND, P2 
XND , TW 
XND,  F 
XND, FBL 
XND ,RHOV2 
XND , SF 
XND , DW 
XND, FREX 
XND, FBLREX 
XND, FBT 
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1074  WRITE  (25,5000)  XND.FBTREX 

1075  RHOV1  =  RHOV2 

1076  RHOl  =  RH02 

1077  W1  =  W2 

1078  MIS  =  M2S 

1079  Ml  =  M2 

1080  U1  =  U2 

1081  PI  =  P2 

1082  2800  CONTINUE 

1083  I  =  28 

1084  PZ  =  P(I) 

1085  PO  =  PCI+1) 

1086  PT  =  P (1+2) 

1087  XZ  =  X(I) 

1088  XO  =  X(I+1) 

1089  XT  =  X  (1+2) 

1090  DO  3000  J  =  1,15 

1091  C 

1092  C  ************************************* 

1093  C  ********  2.580  <  XL  <  2.650  ********* 

1094  C  *****  0.973202  <  XND  <  0.999607  ***** 

1095  C  ************************************* 

1096  C 

1097  XL  =  XL  +  DX 

1098  CALL  SUBPRES  (PZ ,PO, PT ,XZ , XO, XT ,XL , XND,P2) 

1099  IF  (Ml  .LT.  0.01)  THEN 

1100  CALL  SUBINC 

(FUD,A,TO,PEX ,XND,P2 , DPS, RHOl ,RH0V1 ,U1 ,W1 ,M1 , 

1101  +M1S , DWSUM.T , SFSUM.RH02 .RHOV2 , U2 , W2 ,M2S ,M2 , DW ,REX , 

REW, 

1102  +BETA , TW , SF , F , FBL , FBT , PI , FREX , FBLREX , FBTREX , RHOBAR , 

1103  +UBAR , MFFB , MFFS) 

1104  ELSE 

1105  CALL  SUBCOMP 

(FUD , A , TO , PEX , XND , P2 , DPS , PI , RHOl , RHOV 1 , U1 , W1 , 

1106  +M1 , MIS, DWSUM.T, SFSUM.RH02 ,RHOV2 ,U2 ,W2 ,M2S,M2 ,DW, 

REX,  REW, 

1 1 07  +BETA , TW , SF , F , FBL , FBT , FREX , FBLREX , FBTREX , RHOBAR , UBAR) 


1108 

ENDIF 

1109 

WRITE 

(11,5000) 

XND, M2 

1110 

WRITE 

(12,5000) 

XND, REW 

1111 

WRITE 

(13,5000) 

XND, REX 

1112 

WRITE 

(14,5000) 

XND , P2 

1113 

WRITE 

(15,5000) 

XND.TW 

\  1114 

WRITE 

(16,5000) 

XND,  F 

1115 

WRITE 

(17,5000) 

XND, FBL 

1116 

WRITE 

(18,5000) 

XND.RHOV2 

1117 

WRITE 

(19,5000) 

XND , SF 

1118 

WRITE 

(20,5000) 

XND, DW 

1119 

WRITE 

(22,5000) 

XND, FREX 

1120 

WRITE 

(23,5000) 

XND, FBLREX 

1121 

WRITE 

(24,5000) 

XND, FBT 
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1122 

1123 

1124 

1125 

1126 

1127 

1128 

1129 

1130  3000 

1131 

1132 

1133 

1134 

1135 

1136 

1137 

1138 

1139 

1140 

1141 

1142 

1143 

1144 

1145 

1146  C 

1147  C 

1148  C 

1149  C 

1150  C 

1151  4000 

1152  4010 

1153  4020 

1154  4030 

1155  C 

1156  C 

1157  C 

1158  C 

1159  C 

1160  5000 

1161  5010 


1162 

1163  5011 


1164 

1165  5020 

1166  5030 

1167  5040 

1168  5050 

1169  5060 


WRITE  (25,5000)  XND, FBTREX 

RHOV1  =  RHOV2 

RHOl  =  RH02 

VI  =  W2 

MIS  =  M2S 

Ml  =  M2 

U1  =  U2 

PI  =  P2 

CONTINUE 

PI  =  3 . 14 15926536D0 
D  =  0.5/12.0 

CL  =  (31.8125  -  16 . 15625) /12 . 0 

AC  =  PI *D*  *2/4 . 0 

PF  =  (P(l)  -  P (30) ) * AC 

TEST  =  PF  -  SFSUM 

PEF  =  DABS (100. 0*TEST/PF) 

PEM  =  DABS (100 . 0*DWSUM/WCK) 

WRITE  (9,5090)  PF , TEST , PEF ,WCK , DWSUM , PEM 
SFC  =  PF  -  SFEA 
TWBARC  =  SFC/ (PI*D*CL) 

FBARC  =  2 . 0*GC*TWBARC/ (RHOC*UBARC*  *2) 
FCREXC  =  FBARC*R£XBRC 
WRITE  (9,5100)  FBARC 
WRITE  (9,5110)  FCREXC 

************************* 

*****  READ  FORMAT’S  ***** 
************************* 


FORMAT  (ID 
FORMAT  (F15.12) 

FORMAT  (F6.2) 

FORMAT  (F5.3) 

************************** 

*****  WRITE  FORMAT’S  ***** 
************************** 


FORMAT  (1X,F20.13,5X,F20.13) 

FORMAT 

(2X , F4 . 3 , 2X , F8 . 5 , 2X , F9 .4 , 2X , F9 . 3 , 3X , F6 . 4 , 2X , F8 . 6 , 
3X.F9.5, 

+5X , F8 . 5 , 3X , F8 . 5) 

FORMAT 


(2X ,F4 . 3 ,2X ,F8 . 5 ,2X ,F9 .4 ,2X ,F9 . 3 , 3X ,F6 .4 ,2X ,F8 . 6 , 
3X.F9.5, 


+5X.F8.5) 

FORMAT  (45X , ' TEST  RUN 
FORMAT  (15X , ' Patm 
FORMAT  (15X , ' TO 
FORMAT  (15X , ' Psink 
FORMAT  (3X, 'ENTER  THE 


#  ' ,11,////) 

=  ' ,F7.4,lX,'psia' ,/) 

*  ' ,F6.2,3X,’R' ,/) 

=  ’ ,F7.4,lX,’Psia' ,//) 

POROSITY  FLUX  CORRECTION 
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FACTOR,  X.XXX') 

1170  5061  FORMAT  (3X, 'ENTER  THE  BLOWING  MOMENTUM  FLUX 

FACTOR  X . XXX ' ) 

1171  5062  FORMAT’ <3X, 'ENTER  THE  SUCTION  MOMENTUM  FLUX 

FACTOR  X.XXX') 

1172  5070  FORMAT’  (15X,’ Porosity ' ,/,15x, 'Flux  Factor  = 

' , 2X , F5 . 3 , // , 15X , 

1173  +'Blowing  Momentum' ,/, 15X ,' Flux  Factor  = 

' ,2x,F5 .3.//.15X, 

1174  +'Suctir>n  Momentum'  ,/,  15X  ,' Flux  Factor  = 

' , 2X , F5 . 3 , ///) 

1175  5080  FORMAT 

C21X, 'RADIAL' ,6X, 'AXIAL' , 37X , ' FB (REX) ' ,4X, ’FB(REX) 

1176  +  3X ,  ’ X/L ’ ,2X, 'FRICTION' , 3X REYNOLDS ' , 3X REYNOLDS’ 

5X , ' BETA' , 

1177  +6X , ' MACH’ , 6X , ’ F (REX) ' , 6X , ' F .  DEV.’,4X,'F. 

DEV. ' ,/,10X, 'FACTOR' , 

1178  +5X , 'NUMBER' ,5X, 'NUMBER' , 1 3X ,' NUMBER ’ , 18X ,' LAMINAR' 

2X 

1179  + 1  TURBULENT ’  ,/) 

1180  5090  FORMAT  (///, 15X Pressure  Force’, 5X,’= 

’  , FI  8 . 14 , IX , ’ lbf ’  , 

1181  +//, 15X ,’ Force  Balance' , 6X ,' = 

’ ,F18.14,lX,'lbf' ,//,15X,'%  Error' , 

1182  +12X , ' =  ’ ,F8.4,////,15X,'Mass  Flow  Rat e ',/, 15X , 

1183  +’ Entering  Condenser  = 

' , FI 8. 14, IX,’ lbm/sec ' , // , 15X , 'Mass  Balance’ , 

1184  +7X , ' =  ' .F18.14.1X, 'lbm/sec' ,//,15X,’% 

Error’ ,12X, '=  ’ ,F8.4) 

1185  5100  FORMAT  (//, 15X ,' Average 

Condenser’ ,/, 15X ,* Friction  Factor' ,4X, 

1186  +'=  ’ , F18 . 14) 

1187  5110  FORMAT  (//, 15X , 'Average  Condenser ',/, 15X f (Rex) 

product ’  , 5X  , 

1188  +'=  ’ ,F18 . 14) 

1189  C 

1190  CLOSE  (9) 

1191  CLOSE  (10) 

1192  CLOSE  (11) 

1193  CLOSE  (12) 

1194  CLOSE  (13) 

1195  CLOSE  (14) 

1196  CLOSE  (15) 

1197  CLOSE  (16) 

1198  CLOSE  (17) 

1199  CLOSE  (18) 

1200  CLOSE  (19) 

1201  CLOSE  (20) 

1202  CLOSE  (21) 

1203  CLOSE  (22) 

1204  CLOSE  (23) 
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1205  CLOSE  (24) 

1206  CLOSE  (25) 

1207  C 

1208  C  ******************************** 

1209  C  ******  END  OF  MAIN  PROGRAM  ***** 

1210  C  ******************************** 

1211  C 

1212  END 


Name 

Type 

A 

REAL *8 

AC 

REAL* 8 

BETA 

REAL* 8 

CL 

REAL *8 

D 

DABS 

REAL *8 

DPS 

REAL *8 

DW 

REAL* 8 

DWSUM 

REAL *8 

DX 

REAL *8 

F 

REAL *8 

FBARC 

REAL*8 

FBL 

REAL* 8 

FBLREX 

REAL *8 

FBT 

REAL* 8 

FBTREX 

REAL *8 

FCREXC 

REAL* 8 

FREX 

REAL*8 

FUD 

REAL* 8 

GC 

REAL *8 

I 

INTEGER *4 

J 

INTEGER*4 

K 

INTEGER* 4 

Ml 

REAL* 8 

MIS 

REAL *8 

M2 

REAL* 8 

M2S 

REAL *8 

MFFB 

REAL *8 

MFFS 

REAL *8 

N 

INTEGER *4 

P 

REAL *8 

PI 

REAL *8 

P2 

REAL *8 

PATM 

REAL* 8 

PEF 

REAL *8 

PEM 

REAL *8 

PEX 

REAL *8 

PF 

REAL* 8 

PI 

REAL *8 

PL 

REAL* 8 

PND 

REAL *8 

Offset  P  Class 

1034 

1458 

1298 

1450 

1442 

INTRINSIC 

1202 

1274 

1130 

1026 

1322 

1514 

1330 

1354 

1338 

1362 

1522 

1346 

1002 

1114 

998 

1386 

1390 

1066 

1058 

1266 

1258 

1010 

1018 

970 

482 

1082 

1090 

974 

1482 

1490 

1122 

1466 

1434 

962 

722 
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PO 

REAL* 8 

1154 

PS  INK 

REAL *8 

990 

PT 

REAL* 8 

1162 

PZ 

REAL *8 

1146 

R 

REAL *8 

1106 

REW 

REAL *8 

1290 

REX 

REAL* 8 

1282 

REXBRC 

REAL *8 

1402 

RHOl 

REAL* 8 

1210 

RH02 

REAL *8 

1226 

RHOBAR 

REAL *8 

1370 

RHOC 

REAL* 8 

1418 

RHOV1 

REAL *8 

1218 

RHOV2 

REAL *8 

1234 

SF 

REAL *8 

1314 

SFC 

REAL *8 

1498 

SFEA 

REAL *8 

1410 

SFSUM 

REAL *8 

1138 

T 

REAL *8 

1098 

TO 

REAL*8 

982 

TEST 

REAL *8 

1474 

TW 

REAL *8 

1306 

TWBARC 

REAL *8 

1506 

U1 

REAL *8 

1074 

U2 

REAL* 8 

1242 

UBAR 

REAL *8 

1378 

UBARC 

REAL *8 

1394 

W1 

REAL *8 

1050 

W2 

REAL*  8 

1250 

WCK 

REAL *8 

1426 

X 

REAL* 8 

2 

XL 

REAL *8 

1042 

XND 

REAL* 8 

1194 

XNDD 

REAL *8 

242 

XO 

REAL *8 

1178 

XT 

REAL *8 

1186 

xz 

REAL*8 

1170 

1213  C 

1214  C 

**************************************************** 

1215  C 

INTERPOLATION 

!  SUBROUTINE  TO  DETERMINE  LOCAL  PRESSURE 

1216  C 

1217  C 

**************************************************** 

1218 
1219  C 

SUBROUTINE 

:  SUBPRES  (PZ , PO, PT  , X7. , XO, XT  ,XL  , XND , P2) 

1220  C 

**  GIVEN  THREE  PRESSURES,  PZ,  PO,  &  PT,  AT  THREE 
LOCATIONS, 

1221  C 

XZ,  XO,  &  XT,  THIS  SUBROUTINE  INTERPOLATES 

TO  FIND  THE 

1222  C 

PRESSURE, 

<  XT .  XL 

P2 ,  WITH  THE  REQUIREMENT  XO  <  XL 

IS 
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REFERRED  TO  AS  “STATION  2“  ** 


r 


1223  C 

1224  C 

1225  IMPLICIT  REAL *8  (A-H.O-Z) 

1226  PL  =  31 .8125/12.0 

1227  CC1  =  PZ/((XZ  -  X0)*(XZ  -  XT)) 

1228  CC2  =  PO/((XO  -  XZ)*(XO  -  XT)) 

1229  CC3  =  PT/ ( (XT  -  XZ) * (XT  -  XO)  ) 

1230  P2  =  CC1*  (XL  -  XO)  *  (XL  -  XT)  +  CC2*(XL  -  XZ)  *  (XL 
-  XT)  +  CC3*  (XL 

1231  +-  XZ)  *  (XL  -XO) 


1232 

1233 

1234 

XND  = 

RETURN 

END 

XL/PL 

Name 

Type 

Offset 

P  Class 

CC1 

REAL *8 

3076 

CC2 

REAL *8 

3084 

CC3 

REAL *8 

3092 

P2 

REAL *8 

32 

* 

PL 

REAL *8 

3068 

PO 

REAL* 8 

4 

* 

PT 

REAL *8 

8 

* 

PZ 

REAL *8 

0 

* 

XL 

REAL *8 

24 

* 

XND 

REAL *8 

28 

* 

XO 

REAL *8 

16 

* 

XT 

REAL*8 

20 

* 

XZ 

REAL *8 

12 

* 

1235  C 

1236  C  ******************************************* 

1237  C  *****  INCOMPRESSIBLE  MODEL  SUBROUTINE  ***** 

1238  C  ******************************************* 

1239  C 

1240  SUBROUTINE  SUBINC 

(FUD , A , TO , PEX , XND , P2 , DPS , RHOl .RHOV1 ,U1 ,W1 ,M1 , 

1241  +M1S , DWSUM , T , SFSUM . RH02 , RHOV2 , U2 , W2 , M2  S , M2 , DW , REX , 

REW , BETA , T 
W.SF, 

1242 

1242  +F, FBL , FBT , PI , FREX , FBLREX , FBTREX ,RHOBAR , UBAR , MFFB , 

MFFS) 

1243  C 

1244  C  **  THIS  SUBROUTINE  CALCULATES  FLOW  PROPETIES  AT 

1245  C  STATION  2  BASED  ON  AN  INCOMPRESSIBLE  MODEL  ** 

1246  C 

1247  IMPLICIT  REAL *8  (A-H,M,0-Z) 

1248  PI  =  3.14159265  36D0 

1249  D  =  0.5/12.0 

1250  DX  =  0.005 
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1251 

1252 

1253 

1254 

1255 

1256 

1257 

1258 

1259 

1260 
1261 
1262 

1263 

1264 

1265 

1266 

1267 

1268 

1269 

1270 

1271 

1272 

1273 

1274 

1275 

1276 

1277 

1278 

1279 

1280 
1281 
1282 

1283 

1284 

1285 

1286 

1287 

1288 

1289 

1290 

1291 

1292 

1293 

1294 

1295 

1296 

1297 

1298 

1299 

1300 

1301 


AP  =  PI*D*DX 

AC  =  PI*D**2/4 . 0 

GC  =  32.174 

R  =  53.335 

GAM  =1.4 

GAM1  =0.2 

PL  =  31.8125/12.0 

CL  =  (31.8125/2.0  -  0.25)/12.0 

DPS  =  P2*  *2  -  PEX*  *2 

ADPS  =  DABS (DPS) 

RHOV2  =  FUD*GC*ADPS/A 

IF  (XND  .LT.  0.5)  RHOV2  =  -RHOV2 

RHOVAV  =  (RHOV1  +  RHOV2)/2.0 

DW  =  -RHOVAV *AP 

DWSUM  =  DWSUM  +  DW 

W2  =  W1  +  DW 

WBAR  =  (W1  +  W2)/2.0 

RH02  =  P2/(R*T) 

RHOBAR  =  (RHOl  +  RH02)/2.0 

C2  =  GC*GAM*P2/RH02 

U2  =  RH01*U1/RH02  +  DW/(RH02*AC) 

UBAR  =  (U1  +  U2)/2.0 
M2S  =  U2*  *2/C2 
MSBAR  =  (MIS  +  M2S)/2.0 
M2  =  DSQRT (M2S) 

CONI  =  1  +  GAM1*MSBAR 
T  =  TO/CONI 

RMU  =  2.27E-08*DSQRT(T**3)/(T  +  198.6) 
REX  =  4 . Q*WBAR/ (PI *D*RMU*GC) 

REW  =  RHOVAV *D/ (RMU*GC) 

BETA  =  DABS (REW/REX) 

MFF1  =  MFFB 
MFF2  =  MFFB 

IF  (XND  .GT.  0.508)  MFF2  =  MFFS 

IF  (XND  .GT.  0.510)  MFF1  =  MFFS 

DMOM  =  W2*U2*MFF2/GC  -  W1*U1*MFF1/GC 

TW  =  ((PI  -  P2)*AC  -  DMOM) /AP 

SF  =  TW*AP 

SFSUM  =  SFSUM  +  SF 

F  =  2 . 0*GC*TW/ (RHOBAR*UBAR**2) 

AREX  =  DABS (REX) 

FBL  =  16.0/REX* (1 .2337  - 

0 . 2337*DEXP(0 . 0363*REW) ) *DEXP (1 . 2*MSBAR) 

FSTAR  =  0 . 079/AREX**0 .25 

FBT  =  FSTAR* (1.0  +  17 .5*AREX**0 .25* BETA) 

FREX  =  F*REX 

FBLREX  =  FBL* REX 

FBTREX  =  FBT*REX 

IF  (XND  .LT.  0.491)  THEN 

WRITE  (27,6000)  F,FBL 

END  IF 

6000  FORMAT  (IX ,F13 .9 ,2X , F13 .9) 
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1302 

1303 


RETURN 

END 


Name 

Type 

Offset 

P 

Class 

A 

REAL*8 

4 

♦ 

AC 

REAL*  8 

3132 

ADPS 

REAL* 8 

3188 

AP 

REAL*8 

3124 

AREX 

REAL *8 

3268 

BETA 

REAL*8 

100 

* 

C2 

REAL *8 

3212 

CL 

REAL* 8 

3180 

CONI 

REAL *8 

3228 

D 

REAL* 8 

3108 

DABS 

INTRINSIC 

DEXP 

INTRINSIC 

DMOM 

REAL *8 

3260 

DPS 

REAL* 8 

24 

* 

DSQRT 

INTRINSIC 

DW 

REAL *8 

88 

* 

DWSUM 

REAL *8 

52 

* 

DX 

REAL *8 

3116 

F 

REAL *8 

112 

* 

FBL 

REAL*  8 

116 

* 

FBLREX 

REAL *8 

132 

* 

FBT 

REAL *8 

120 

* 

FBTREX 

REAL *8 

136 

* 

FREX 

REAL*  8 

128 

♦ 

FSTAR 

REAL*  8 

3276 

FUD 

REAL*8 

0 

* 

GAM 

REAL*  8 

3156 

GAM1 

REAL* 8 

3164 

GC 

REAL *8 

3140 

Ml 

REAL* 8 

44 

* 

MIS 

REAL *8 

48 

* 

M2 

REAL *8 

84 

* 

M2S 

REAL *8 

80 

4c 

MFFl 

REAL *8 

3244 

MFF2 

REAL*8 

3252 

MFFB 

REAL *8 

148 

* 

MFFS 

REAL*8 

152 

4c 

MS  BAR 

REAL *8 

3220 

PI 

REAL* 8 

124 

4t 

P2 

REAL *8 

20 

4c 

PEX 

REAL*8 

12 

* 

PI 

REAL* 8 

3100 

PL 

REAL *8 

3172 

R 

REAL* 8 

3148 

REW 

REAL *8 

96 

4c 

REX 

REAL *8 

92 

* 

RHOl 

REAL *8 

28 

4c 
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RH02 

REAL* 8 

64 

* 

RHOBAR 

REAL *8 

140 

* 

RHOV1 

REAL *8 

32 

♦ 

RHOV  2 

REAL *8 

68 

* 

RHOVAV 

REAL* 8 

3196 

RMU 

REAL* 8 

3236 

SF 

REAL* 8 

108 

* 

SFSUM 

REAL *8 

60 

♦ 

T 

REAL* 8 

56 

* 

TO 

REAL *8 

8 

* 

TW 

REAL *8 

104 

* 

U1 

REAL *8 

36 

* 

U2 

REAL* 8 

72 

* 

UBAR 

REAL *8 

144 

* 

W1 

REAL*  8 

40 

* 

W2 

REAL *8 

76 

* 

WBAR 

REAL* 8 

3204 

XND 

REAL *8 

16 

* 

1304  C 

1305  C 

1306  C 

1307  C 

1308  C 

1309 

1310 

1311 

1312  C 

1313  C 

1314  C 

1315  C 

1316 

1317 

1318 

1319 

1320 

1321 

1322 

1323 

1324 

1325 

1326 

1327 

1328 

1329 

1330 

1331 

1332 

1333 


*****  COMPRESSIBLE  MODEL  SUBROUTINE  ***** 


SUBROUTINE  SUBCOMP 

(FUD , A , TO , PEX , XND , P2 , DPS , PI , RHOl , RHOV 1 , U1 , W1 , 

+M1 ,M1  S , DWSUM ,T , SFSUM , RH02 , RHOV2 , U2 , W2 , M2  S , M2 , DW , 
REX , REW , BET , A , TW , 

+SF,F,FBL, FBT , FREX , FBLREX , FBTREX ,RHOBAR , UBAR) 

**  THIS  SUBROUTINE  CALCULATES  FLOW  PROPETIES  AT 
STATION  2  BASED  ON  A  COMPRESSIBLE  MODEL  ** 

IMPLICIT  REAL *8  (A-H,M,0-Z) 

PI  =  3 . 14 15  926536D0 

D  =  0.5/12.0 

DX  =  0.005 

AP  =  PI *D*DX 

AC  =  PI*D**2/4 . 0 

GC  =  32.174 

R  =  53.335 

GAM  =1.4 

GAM1  =0.2 

PL  =  31.8125/12.0 

CL  =  (31.8125/2.0  -  0.25)/12.0 

DPS  =  P2 * *2  -  PEX* *2 

ADPS  =  DABS (DPS) 

RHOV 2  =  FUD*GC*ADPS/A 
IF  (XND  .LT.  0.5)  RHOV 2  =  -RHOV2 
RHOVAV  =  (RHOV1  +  RHOV2)/2.0 
DW  =  -RHOVAV *AP 
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1334  DWSUM  =  DWSUM  +  DW 

1335  W2  =  W1  +  DW 

1336  WBAR  =  (W1  +  W2>/2.0 

1337  CM  =  MIS* (1.0  + 

GAM1*M1S) * (CW2/W1) * *2) / ( (P2/P1) * *2) 

1338  M2S  =  (-1.0  +  DSQRT(1.0  + 

4 . 0*GAM1*CM) ) / (2 . 0*GAM1) 

1339  MSBAR  =  (MIS  +  M2S)/2.0 

134  0  M2  =  DSQRTCM2S) 

1341  CONI  =  1  +  GAM1*MSBAR 

1342  T  =  TO/CONI 

134  3  RH02  =  P2/(R*T) 

1344  RHOBAR  =  (RHOl  +  RH02)/2.0 

1345  C2  =  GC*GAM*P2/RH02 

1346  U2  =  M2*DSQRT(C2) 

1347  UBAR  =  (U1  +  U2>/2.0 

1348  RMU  =  2 . 27E-08*DSQRT (T**  3) /  (T  +  198.6) 

1349  REX  =  4 . 0*WBAR/ (PI*D*RMU*GC) 

1350  REW  =  RHOVAV*D/ (RMU*GC) 

1351  BETA  =  DABS (REW/REX) 

1352  CFP  =  (1.0  +  GAM*MSBAR) / (GAM*MSBAR) 

1353  F  =  0.5* (D/DX) * (DLOG (M1S/M2S)  -  CFP*DLOG (P2/P1) ) 

1354  TW  =  0 . 5*F*RHOBAR*UBAR*  *2/GC 

1355  SF  =  TW*AP 

1356  SFSUM  =  SFSUM  +  SF 

1357  AREX  =  DABS (REX) 

1358  FBL  =  16.0/REX* (1 .2337  - 

0 . 2  337*DEXP (0 . 0363*REW) ) *DEXP (1 . 2  *MSBAR) 

1359  FSTAR  *  0 . 079/ AREX* * 0 . 25 

1360  FBT  =  FSTAR* (1.0  +  17 ,5*AREX**0 .25* BETA) 

1361  FREX  =  F*REX 

1362  FBLREX  =  FBL*REX 

1363  FBTREX  =  FBT* REX 

1364  IF  (XND  .LT.  0.491)  THEN 

1365  WRITE  (27,6000)  F.FBL 

1366  ENDIF 

1367  6000  FORMAT  (IX , FI 3 . 9 , 2X , FI 3 . 9) 

1368  RETURN 

1369  END 


Name  Type 

A  REAL* 8 
AC  REAL *8 
ADPS  REAL *8 
AP  REAL *8 
AREX  REAL* 8 
BETA  REAL *8 
C2  REAL* 8 
CFP  REAL *8 
CL  REAL *8 
CM  REAL *8 


Offset  P  Class 

4  * 

3336 
3392 
3328 
34  64 
104  * 

3440 

3456 

3384 

3416 
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CONI 

REAL *8 

34  32 

D 

DABS 

DEXP 

DLOG 

REAL *8 

3312 

INTRINSIC 

INTRINSIC 

INTRINSIC 

DPS 

DSQRT 

REAL  *  8 

24  * 

INTRINSIC 

DW 

REAL *8 

92  * 

DWSUM 

REAL *8 

56  * 

DX 

REAL *8 

3320 

F 

REAL* 8 

116  * 

FBL 

REAL *8 

120  * 

FBLREX 

REAL*8 

132  * 

FBT 

REAL *8 

124  * 

FBTREX 

REAL* 8 

136  * 

FREX 

REAL* 8 

128  * 

FSTAR 

REAL* 8 

3472 

FUD 

REAL*8 

0  * 

GAM 

REAL *8 

3360 

GAM1 

REAL *8 

3368 

GC 

REAL*  8 

3344 

Ml 

REAL *8 

48  * 

MIS 

REAL  *8 

52  * 

M2 

REAL *8 

88  * 

M2S 

REAL  *8 

84  * 

MSBAR 

REAL *8 

3424 

PI 

REAL *8 

28  * 

F2 

REAL *8 

20  * 

PEX 

REAL* 8 

12  * 

PI 

REAL *8 

3304 

PL 

REAL *8 

3376 

R 

REAL *8 

3352 

REW 

REAL* 8 

100  * 

REX 

REAL *8 

96  * 

RHOl 

REAL *8 

32  * 

RH02 

REAL *8 

68  * 

RHOBAR  REAL *8 

140  * 

RHOV1 

REAL *8 

36  * 

RHOV2 

REAL* 8 

72  * 

RHOVAV  REAL *8 

3400 

RMU 

REAL* 8 

3448 

SF 

REAL *8 

112  * 

SFSUM 

REAL* 8 

64  * 

T 

REAL *8 

60  * 

TO 

REAL* 8 

8  * 

TW 

REAL *8 

108  * 

U1 

REAL*8 

40  * 

U2 

REAL *8 

76  * 

UBAR 

REAL* 8 

144  * 

W1 

REAL *8 

44  * 

W2 

REAL* 8 

80  * 

WBAR 

REAL *8 

3408 
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Name 


Type 


Size 


Class 


ONED 

SUBCOM 

SUBINC 

SUBPRE 


PROGRAM 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 


No  Errors  Detected 
1369  Source  Lines 


Pass  One 


Appendix  C 
TEST  RUN  #  1 


Patm  =  14.2827  psia 

TO  =  526.50  R 
Psink  =  14.2444  psia 

Porosity 

Flux  Factor  =  1.000 

Blowing  Momentum 
Flux  Factor  =  1.290 


Suction  Momentum 
Flux  Factor  =  1.470 


RADIAL 

X/L  FRICTION  REYNOLDS 

FACTOR 

NUMBER 

.100 

.  18663 

-3.4595 

.200 

.09210 

-3.4621 

.300 

.06123 

-3.4665 

.400 

.04737 

-3.4726 

.600 

.03083 

3.4723 

.700 

.05958 

3.4734 

.800 

.08847 

3.4740 

.900 

. 12791 

3.4745 

Pressure 

Force 

= 

Force  Balance 

= 

%  Error 

= 

Mass  Flow  Rate 

Entering 

Condenser  = 

Mass  Balance 

= 

%  Error 

= 

Average 

Condenser 

Friction 

Factor 

= 

Average 

Condenser 

f  (Rex)  product 

*  13 

AXIAL 

REYNOLDS 

BETA 

MACH 

NUMBER 

NUMBER 

87.995 

.0393 

.000312 

176.035 

.0197 

.000625 

264.163 

.0131 

.000938 

352  .423 

.0099 

.001251 

354.203 

.0098 

.001252 

265.850 

.0131 

.000938 

177.478 

.0196 

.000624 

89.092 

.0390 

.000310 

00001883160918  lbf 
00000000005909  lbf 
0003 

00017256225293  lbm/sec 
00000005382849  lbm/sec 
0312 

06447082216418 

99607226000200 


F  (REX) 


16.42282 
16.21332 
16.17488 
16.69325 
10.91849 
15.84048 
15 .70187 
11.39598 
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TEST  RUN  #  2 


Patm  =  14.2484  psia 

TO  =  526.00  R 
Psink  =  14.2281  psia 

Porosity 


Flux 

Factor  = 

1.000 

Blowing  Momentum 

Flux 

Factor  = 

1.305 

Suction  Momentum 

Flux 

Factor  = 

1.370 

RADIAL 

X/L 

FRICTION 

REYNOLDS 

FACTOR 

NUMBER 

.  100 

.34820 

-1.8347 

.200 

.17400 

-1.8357 

.300 

.11622 

-1.8375 

.400 

.08698 

-1.8399 

.600 

.08404 

1.8414 

.700 

.11145 

1.8409 

.800 

.16664 

1.8405 

.900 

.33160 

1.8403 

Pressure  Force  = 

Force  Balance  = 

%  Error  = 

Mass  Flow  Rate 
Entering  Condenser  = 

Mass  Balance  = 

%  Error  = 

Average  Condenser 
Friction  Factor  = 


AXIAL 

REYNOLDS 

BETA 

MACH 

NUMBER 

NUMBER 

46.669 

.0393 

.000166 

93.355 

.0197 

.000332 

140.076 

.0131 

.000498 

186.850 

.0098 

.000664 

187.725 

.0098 

.000664 

140.887 

.0131 

.000497 

94.061 

.0196 

.000331 

47.242 

.0390 

.000165 

00001122276190  lbf 
00000000000089  lbf 
0000 

00009140706303  lbm/sec 
00000000598701  lbm/sec 
0065 

13496525522116 


Average  Condenser 

f (Rex)  product  =  15.53163701557668 


F (REX) 


16.24985 
16.24329 
16.27910 
16.25176 
15.77694 
15 .70123 
15.67422 
15.66552 


C-2 


TEST  RUN  #  3 


Patm  =  14.3022  psia 

TO  =  526.00  R 
Psink  =  14.2308  psia 


Porosity 

Flux  Factor  =  1.000 

Blowing  Momentum 
Flux  Factor  =  1.270 

Suction  Momentum 
Flux  Factor  »  1.600 


RADIAL 

AXIAL 

X/L 

FRICTION 

REYNOLDS 

REYNOLDS 

BETA 

MACH 

F  (REX) 

FACTOR 

NUMBER 

NUMBER 

NUMBER 

.100 

.10239 

-6  .4532 

164 . 129 

.0393 

.000582 

16.80494 

.200 

.05118 

-6.4602 

328.378 

.0197 

.001164 

16.80719 

.300 

.03387 

-6.4720 

492.866 

.0131 

.001747 

16.69245 

.400 

.02542 

-6.4885 

657.713 

.0099 

.002332 

16.72037 

.600 

.02416 

6.4716 

661.560 

.0098 

.002334 

15.98275 

.700 

.02319 

6.4815 

496.794 

.0130 

.001750 

11.52303 

.800 

.04039 

6.4884 

331.809 

.0196 

.001165 

13.40302 

.900 

.06137 

6.4937 

166.669 

.0390 

.000580 

10.22789 

Pressure  Force  = 

Force  Balance  = 

%  Error  = 

Mass  Flow  Rate 
Entering  Condenser  * 

Mass  Balance  = 


.00002995839503  lbf 
.00000000004151  lbf 
.0001 

.00032190622080  lbm/sec 
-.00000002620157  lbm/sec 


%  Error 


.0081 


Average  Condenser 

Friction  Factor  *  .03506096928445 


Average  Condenser 

f (Rex)  product  =  14.20926769146199 
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TEST  RUN  #  4 


Patm  =  14.2729  psia 

TO  =  526.00  R 
Psink  =  14.1330  psia 


Porosity 

Flux  Factor  =  1.000 

Blowing  Momentum 
Flux  Factor  =  1.245 

Suction  Momentum 
Flux  Factor  =  1.600 


RADIAL 

AXIAL 

X/L 

FRICTION 

REYNOLDS 

REYNOLDS 

BETA 

MACH 

F (REX) 

FACTOR 

NUMBER 

NUMBER 

NUMBER 

.100 

.05462 

-12.5627 

319.471 

.0393 

.001138 

17.45062 

.200 

.02702 

-12.5848 

639.317 

.0197 

.002277 

17.27199 

.300 

.01815 

-12.6219 

959.915 

.0131 

.003418 

17.41907 

.400 

.01360 

-12.6740 

1281.648 

.0099 

.004564 

17.43204 

.600 

.01261 

12.5920 

1290.864 

.0098 

.004575 

16.27655 

.700 

.01401 

12.6353 

969.950 

.0130 

.003432 

13.58974 

.800 

.02095 

12.6676 

648.078 

.0195 

.002286 

13.57851 

.900 

.12140 

12.6756 

325.702 

.0389 

.001138 

39.54022 

Pressure  Force  = 

Force  Balance  = 

%  Error  = 

Mass  Flow  Rate 
Entering  Condenser  = 

Mass  Balance  = 

%  Error  = 

Average  Condenser 
Friction  Factor  = 

Average  Condenser 
f (Rex)  product  = 


.00005293412287  lbf 
.00000000069008  lbf 
.0013 

.00062766683506  lbm/sec 
.00000004480944  lbm/sec 
.0071 

.02606307347730 


20.59576260070104 
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